0JC  1SBEJS 

wwasiiaaif 


THE  UNIVERSITY  OF  ALBERTA 


RELEASE  FORM 

NAME  OF  AUTHOR  MOHAMMED  ZAHEERUDDIN 

TITLE  OF  THESIS  COMPUTER  SIMULATION  OF  DIRECT  GAIN 

PASSIVE  HOUSE 

DEGREE  FOR  WHICH  THESIS  WAS  PRESENTED  MASTER  OF  SCIENCE 
YEAR  THIS  DEGREE  GRANTED  FALL  1980 

Permission  is  hereby  granted  to  THE  UNIVERSITY  OF 
ALBERTA  LIBRARY  to  reproduce  single  copies  of  this 
thesis  and  to  lend  or  sell  such  copies  for  private, 
scholarly  or  scientific  research  purposes  only. 

The  author  reserves  other  publication  rights,  and 
neither  the  thesis  nor  extensive  extracts  from  it  may 
be  printed  or  otherwise  reproduced  without  the  author's 
written  permission. 


THE  UNIVERSITY  OF  ALBERTA 


COMPUTER  SIMULATION  OF  DIRECT  GAIN  PASSIVE  HOUSE 


by 

MOHAMMED  ZAHEERUDDIN 


A  THESIS 

SUBMITTED  TO  THE  FACULTY  OF  GRADUATE  STUDIES  AND  RESEARCH 
IN  PARTIAL  FULFILMENT  OF  THE  REQUIREMENTS  FOR  THE  DEGREE 

OF  MASTER  OF  SCIENCE 

DEPARTMENT  OF  MECHANICAL  ENGINEERING 


EDMONTON,  ALBERTA 


FALL  1980 


THE  UNIVERSITY  OF  ALBERTA 


FACULTY  OF  GRADUATE  STUDIES  AND  RESEARCH 

The  undersigned  certify  that  they  have  read,  and 
recommend  to  the  Faculty  of  Graduate  Studies  and  Research, 
for  acceptance,  a  thesis  entitled  COMPUTER  SIMULATION  OF 
DIRECT  GAIN  PASSIVE  HOUSE  submitted  by  MOHAMMED  ZAHEERUDDIN 
in  partial  fulfilment  of  the  requirements  for  the  degree  of 
MASTER  OF  SCIENCE. 


ABSTRACT 


An  interactive  computer  program  to  predict  the  thermal 
performance  of  any  real  house  was  developed.  An  improved 
method  of  calculating  the  solar  heat  gain  through  single, 
double,  or  triple  glazing  windows  was  incorporated  in  the 
model.  The  winter  solar  gains  can  be  maximized  by  sloping 
the  windows.  Two  different  techniques  of  reducing  summer 
overheating  were  studied.  Either  heat  absorbing  glass 
windows  or  archi tectura 1  projections  such  as  overhangs  may 
be  used. 

Two  important  considerations  in  modelling  the  behaviour 
of  a  house  with  thermal  mass  wall  are  :  the  optimum 
thickness  of  the  thermal  mass  wall  and  the  effect  of  the 
number  of  nodes  on  the  overall  accuracy  of  the  annual 
heating  load  calculations.  It  was  found  that  thermal  mass  in 
the  form  of  5  inch  thick  masonry  wall  exposed  on  both  sides 
to  room  air  or  masonry  in  the  form  of  2.5  inch  thick 
surfacings  on  the  interior  walls  was  most  efficient.  In 
modelling  the  dynamics  of  the  thermal  mass  wall  by  finite 
difference  calculations,  division  of  mass  wall  into  1/2  inch 
thick  layers  ensures  the  accuracy  of  the  annual  heating  load 
calculations  within  0.3%  of  the  calculations  carried  out 
with  twice  the  number  of  nodes  in  the  same  thickness  of  mass 
wall. 

Sizing  of  the  windows  depending  on  the  time  constant  of 
the  house  can  be  studied  from  the  general  design  curves.  The 
applicability  of  the  design  curves  irrespective  of  the 


TV 


weather  data  and  the  internal  heat  load  schedule  was 
demonstrated.  Two  important  parameters,  the  SIG  factor  and 
the  AHLR  factor,  along  with  the  time  constant,  fully 
describe  the  thermal  design. 

The  SIG  factor  gives  a  quantitative  asessment  of  the 
heat  potential  avialable  from  solar  and  the  fluctuating 
component  of  internal  load  that  occur  during  the  day.  The 
AHLR  factor  is  an  output  parameter  which  describes  how  much 
auxiliary  heat  is  required  in  the  house. 

The  economics  of  the  passive  design  features 
(insulation  level  of  the  house,  window  shutter  system  and 
the  thermal  capacity  in  the  house)  were  examined  by 
employing  the  life  cycle  cost  analysis. 
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1.  INTRODUCTION 


Solar  passive  heating  is  one  of  the  most  promising 
alternatives  avialable  to  reduce  house  heating  requirements. 
In  the  context  of  building  heating  a  definition  of  passive 
solar  heating  system  was  formulated  by  the  Department  of 
Housing  and  Urban  Development  (HUD)  which  states: 

A  passive  solar  energy  system  is  one  which  uses  the 
building  structure  as  a  collector,  storage,  and 
transfer  mechanism  with  a  minimum  amount  of 
mechanical  equipment. 

Many  attempts  have  been  made  to  define  passive  solar  heating 
system  but  none  is  completely  satisfactory  because  of  the 
wide  variety  of  the  forms  which  such  systems  can  take.  The 
type  of  design  that  was  studied  here  is  usually  called 
direct  gain.  That  is  solar  heating  occurs  directly  through 
the  windows  of  the  house  and  no  separate  form  of  solar 
col  lector  exists. 

It  will  be  useful  to  evolve  a  concept  of  an  ideal  solar 
passive  heating  system.  An  ideal  solar  passive  system  can  be 
looked  upon  as  one  that  requires  no  auxiliary  energy  to 
maintain  inside  comfort  conditions.  Based  on  such  an  ideal 
concept,  it  can  be  said  that  a  building  structure  should  be 
so  designed  that  it  controls  the  amount  of  solar  energy 
reaching  the  interior  masses  of  the  building,  directly  as 
sunlight,  and  transacting  all  heat  exchanges  and  storage  at 
the  building  interior  masses.  In  effect  such  a  building 
would  behave  like  a  temperate  solar  collector.  With  this 
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approach  part  of  the  external  structure  must  be  composed  of 
a  transparent  insulation  (affording  visibility),  and  at  the 
same  time  the  total  external  structure  should  act  as  a  sun 
light  modulator.  Although  the  ideal  solar  passive  system  may 
not  be  achievable  in  practice  one  can  approach  this  ideal  in 
a  building  if  the  thermal  components  of  the  building  have 
been  designed  as  a  single  entity  right  from  the  outset. 

Although  the  solar  passive  system  has  advantages  of 
being  simple  and  economic,  the  calculations  involved  in  the 
design  are  relatively  sophisticated.  This  stems  from  the 
fact  that  a  complex  pattern  of  heat  exchanges  take  place  in 
and  around  any  real  enclosure. 

Today's  architects  and  designers  are  really  faced  with 
the  problems  of  determining  parameters  like:  Resistance 
values  of  enclosure  elements,  amount  of  heat  capacity  in  the 
structure,  the  area  of  window  as  well  as  its  orientation, 
and  the  operational  characteristics  of  window  shutters,  the 
design  of  overhang  etc.  and  eventually  they  are  faced  with 
the  problem  of  determining  the  annual  auxiliary  heat 
requirements  for  their  designs. 

It  is  the  purpose  of  the  present  work  to  provide  the 
designers  and  researchers  with  an  interactive  computer 
program  that  can  accept  every  possible  design  variable(s) 
and  can  model  their  design  by  simulating  the  dynamic 
response  of  the  house  to  temporal  variations  in  radiation 
fluxes  and  ambient  air  temperatures . 

Any  thermal  model  is  necessarily  only  an  approximation 
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to  the  complex  heat  exchanges  taking  place  in  and  around  any 
real  enclosure.  Then  the  question  arises  as  to  how 
accurately  the  thermal  model  predicts  the  performance  of  any 
real  enclosure. 

The  above  subject  of  calculation  of  energy  consumption 
in  residential  buildings  using  simulation  techniques  has 
been  studied  in  recent  years  by  Konard  and  Larsen  (1), 
Zulfikar  et.al  (2)  ,  Davies  (3),  Balcomb  et  a  1 . ( 4 ) ,  Mazria 
and  Wessling  (5),  DOE  -  2  "Energy  Analysis  Program"  (6), 
etc.  On  the  other  hand  the  design  and  performance  of  direct 
gain  passive  houses  using  static  calculations  has  been 
studied  by  many  investigators,  Dumont,  Besant  and  Green  (7), 
Jones  and  Tymura(8),  Balcomb  and  MacFarland  (9)  etc.  The 
performance  of  passive  test  units  has  been  reported  by 
Wheeling  et  al.  (10),  Dumont  and  Besant  et  a  1.(11),  Chen 
et.al.  (12),  Palmiter  et  al.  (13),  Balcomb  et  al.(4),  and 
Jones  (14)  etc . 

However,  the  above  simulation  models  were  either  too 
specific  or  too  complex  requiring  prohibitively  high 
computer  (CPU)  time  to  study  the  design  of  passively  heated 
houses . 

Gilpin, (15)  developed  a  computer  model  of  a  house 
employing  a  south  facing  window  and  a  mass  wall.  An 
essential  feature  of  this  model  is  that  it  is  a  dynamic 
model . 

The  present  study  is  an  extension  of  Gilpin's  thermal 
model  of  passive  house.  An  improved  method  of  calculating 
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the  solar  heat  gain  was  developed.  Greater  emphasis  was 
placed  on  two  important  elements  of  the  passive  house 
namely:  the  windows  and  the  mass  walls. 

The  logical  question  then  arises  as  to  what  is  the  best 
combination  of  the  windows  and  the  mass  wall  for  a  given 
floor  area  that  results  in  the  least  auxiliary  heat 
requirement  for  maintaining  the  inside  comfort  conditions. 

Several  alternatives  are  available.  Collector  windows 
could  be  sloped  to  maximize  the  solar  heat  gain.  Triple 
glazing  or  heat  absorbing  glass  windows  could  be  used  to 
minimize  the  summer  overheating.  Overhangs  could  be  provided 
on  the  sloping  windows. 

Two  different  techniques  of  storing  heat  in  the  mass 
wall  can  be  studied.  In  the  first  technique  the  mass  wall 
can  be  placed  inside  the  house  and  the  heat  flow  into  and 
out  of  the  wall  can  be  studied  as  a  result  of  rise  or  fall 
in  the  room  temperature.  In  the  second  technique  the  mass 
wall  could  be  oriented  to  receive  the  direct  solar  heat  gain 
through  the  window  as  in  the  case  of  Trombe  Wall. 

In  modelling  the  behaviour  of  the  mass  walls,  stability 
and  accuracy  of  the  finite  difference  calculations  is 
important.  One  is  often  faced  with  problems  like:  What  is 
the  optimum  thickness  of  the  mass  wall  for  diurnal  energy 
storage?  How  few  nodes  could  be  considered  in  studying  the 
dynamics  of  the  mass  wall  without  impairing  the  accuracy? 

The  present  work  examines  the  above  alternatives  and 
provides  logical  solutions  to  the  possibility  of  minimum 
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energy  for  maximum  comfort  through  solar  passive  design. 

It  was  recognized  that  not  all  the  designers  and 
architects  have  access  to  computing  facility.  As  a  result 
they  cannot  take  full  advantages  offerred  by  the  present 
interactive  computer  program.  To  help  such  designers,  a  set 
of  general  design  curves  have  been  prepared.  These  design 
curves  are  the  useful  means  of  estimating  the  annual  energy 
consumption  of  any  house. 

Alternatively  parameters  like,  window  and  mass  wall 
areas,  for  minimum  auxiliary  heat  requirement,  can  also  be 
estimated  from  the  design  curves.  These  design  curves  could 
be  used  irrespective  of  the  weather  data  and  the  internal 
heat  load  schedule.  This  is  possible  with  an  additional 
table  on  the  ISHG ( Integrated  solar  heat  gain)  factors  for 
various  locations. 

These  ideas  and  the  accompanying  solutions  are 
presented  in  chapters  2  through  5.  The  physical  model  and 
the  governing  equations  are  discussed  in  chapters  2  through 
4.  Chapter  5  presents  the  flow  charts  illustrating  the 
solution  scheme.  Some  results  like  :  the  transmittance  of 
window  glass  as  a  function  of  angle  of  incidence,  the 
optimum  mass  wall  thickness  for  diurnal  energy  storage,  the 
optimum  number  of  nodes  for  modelling  the  dynamics  of  the 
mass  wall,  the  general  design  curves  and  the  economics  of 
the  passive  design  features  are  presented  in  chapter  6.  This 
is  followed  by  the  conclusions  drawn  from  this  study. 


2.  THE  PHYSICAL  MODEL 


2 . 1  Introduct ion 

A  mathematical  model  has  been  developed  to  enable 
predictions  to  be  made  of  the  thermal  behaviour  of  any  real 
house.  An  essential  feature  of  the  model  is  its  flexibility 
in  modelling  the  dynamic  behaviour  of  almost  any  type  of 
design  and  construction.  Some  of  the  flexibilities  offered 
by  the  model  are: 

1.  Provision  to  model  solar  heat  gain  from  any  number  of 
windows  (as  in  the  real  house)  facing  north,  south, 
east ,  or  west . 

2.  Orientation  of  the  house  could  be  varied  to  study  the 
optimum  house  orientation  angle. 

3.  Solar  heat  gain  can  be  maximized  by  sloping  one  or  all 
of  the  windows  as  desired. 

4.  Summer  overheating  can  be  reduced  by  judicious  design  of 
overhang  on  the  windows. 

5.  Another  useful  technique  of  reducing  the  overheating  can 
be  studied  by  varying  the  number  of  glazings  and  their 
KL  product.  A 1 ternat i vely , the  effectiveness  of  heat 
absorbing  glass  in  reducing  the  overheating  can  also  be 
studied . 

6.  The  possibility  of  using  automatic  shutters  which  open 
and  close  according  to  the  thermostatic  control  switch 
exists  in  the  model.  Alternatively,  shutters  which  open 
at  sunrise  and  close  at  sunset  during  winter  months, 
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while  being  left  open  all  through  summer  months  can  be 
examined.  The  performance  of  both  types  of  shutters  can 
be  compared  with  the  case  of  no  shutters. 

7.  The  model  can  investigate  such  possibilities  as 

placement  of  the  thermal  mass  exterior  or  interior  to 
the  structure.  Further  the  possibility  of  modelling  the 
Trombe  type  of  storage  wall  is  also  included. 

A  number  of  such  possibilities  can  be  studied  from  the 
present  thermal  model.  The  theory  and  the  calculations  are 
illustrated  as  and  when  they  appear  in  the  following 
chapters . 


2 . 2  Phys i ca 1  Mode  1 

2.2.1  The  Therma 1  Enclosure 

Figure  2.1  shows  the  schematic  of  one  of  the  several 
house  designs  the  present  model  can  handle.  The  thermal 
enclosure  1  is  comprised  of  elements  such  as  north,  south, 
east  and  west  facing  walls,  above  and  below  grade  walls, 
ceiling  and  basement  floor.  The  exterior  walls  of  the 
enclosure  are  assumed  to  be  of  standard  construction.  The 
possibility  of  modelling  the  heavier  construction  (such  as 
brick  or  concrete  walls)  in  place  of  wood  frame  construction 
can  also  be  analysed. 
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Figure  2.1  Schematic  of  a  house  illustrating  the  design  components 


9 


A  single  overhang  6  is  an  architectural  projection  for 
shading  the  sloping  windows  against  high  altitudes  of  summer 
radiation. 

The  whole  orientation  of  the  building  could  be  varied 
with  respect  to  true  south  as  shown  in  figure  2.2. 

Within  the  thermal  enclosure,  internal  heat  gains  due 
to  lights,  appliances  and  people  were  considered.  Based  on 
experience  and  literature  survey  (2)  a  reasonable  schedule 
of  internal  heat  gains  was  worked  out.  This  is  shown  in 
figure  2.3. 
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Figure  2.2  The  definition  of  the  house  orientation  angle 
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Figure  2.3  Internal  heat  load  schedule 
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2.2.2  The  W i ndow  ^  Shutter  System 

Most  of  the  problems  associated  with  passive  design 
such  as  heat  loss  in  winters  (apart  from  contributions  like 
heat  loss  due  to  conduction  and  infiltration),  and 
overheating  in  summers  are  attributable  to  the  window 
element;  therefore,  its  design  plays  an  important  part  in  a 
direct  gain  passive  system. 

Figure  2.4  shows  a  sloped  triple  glazing  window  facing 
due  south.  The  possibility  of  partial  shading  of  window  can 
be  studied  by  providing  a  solid  overhang  6.  Large  south 
facing  windows  though  desirable  from  the  solar  energy 
collection  point  of  view,  introduce  two  major  problems.  In 
the  first  place,  during  winters,  the  large  windows  could  be 
net  losers  of  heat  in  the  absence  of  night  insulating 
shutters.  Secondly  they  would  admit  large  solar  radiation  in 
summer  months  which  could  result  in  overheating. 

In  order  to  prevent  the  net  loss  of  heat  and  to  control 
the  excessive  heat  gain,  two  systems  were  proposed.  In  the 
first  case  it  was  assumed  that  a  temperature  sensor 
depending  on  the  room  temperature,  would  transmit  signals  to 
a  shutter  mechanism  to  open  or  close  the  shutters,  in 
proportion  to  the  control  signal.  The  problem  of  admitting 
day  light  irrespective  of  net  loss  or  gain  of  heat  was 
incorporated  in  the  model  by  specifying  a  minimum  opening  of 
shutters  during  the  day  light  hours. 
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Figure  2.4  The  geometry  of  a  sloping  window  and  roof 

overhang 
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The  above  mechanism  is  somewhat  sophisticated  and  to 
date  no  such  mechanism  has  been  developed  commercially. 
Another  logical  and  simple  shutter  routine  could  be  thought 
of  as  one  that  opens  at  sun  rise  and  closes  at  sunset  during 
the  winter  months.  The  shutters  could  be  left  open  during 
summer  months.  The  shutter  operating  schedule  is  depicted  in 
figure  2.5.  A  reasonable  assumption  that  goes  with  this 
shutter  routine  is  that  the  people  would  do  something  when 
the  inside  room  temperature  exceeds  the  maximum  value.  This 
can  be  expected  in  summer  months.  It  is  natural  to  believe 
that  people  would  pull  the  drapes,  close  the  Venetian  blinds 
or  open  the  windows  when  they  sense  overheating. 

Hence  provision  was  made  in  the  model  for  comparing  the 
performance  of  the  above  two  shutter  systems  with  that  of  a 
case  with  no  shutters. 


■ 
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Figure  2.5  Manual  shutter  operating  schedule 
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2.2.3  The  Therma 1  Mass 

The  thermal  mass  in  the  structure  was  assumed  to  occur 
in  one  of  the  two  forms.  That  is  the  thermal  mass  could 
consist  either  of  thin  materials  such  as  plaster  board  or 
thick  materials  such  as  masonry  walls.  In  the  thermal  model 
the  thin  plaster  board  was  treated  by  a  lumped-capacity 
approach.  Because  it  is  reasonable  to  assume  the  thin 
plaster  board  to  be  at  uniform  temperature.  Applying  the 
same  analogy  the  thermal  capacity  introduced  by  the  wood 
frame  construction  as  well  as  the  furniture  in  the  house  is 
studied  by  Jones  and  Tymura  (8).  However,  this  has  to  be 
approached  with  caution  because  not  all  the  wood  used  in  the 
house  directly  interacts  with  the  inside  room  air.  The  thick 
masonry  walls  shown  by  the  elements  7  and  8  in  figure  2.1 
have  to  be  modelled  using  a  distributed  -  capacitance 
approach.  This  is  necessary  because  a  diurnal  temperature 
variation  will  occur  through  a  thick  mass  wall.  This 
phenomenon  was  modelled  using  a  finite  difference 
calculation. 

The  masonry  wall  8  could  be  considered  to  exist  in  one 
of  the  two  forms.  Either  as  a  wall  of  thickness  (t)  exposed 
on  both  sides  to  the  room  air,  or  as  masonry  surfacings  of 
thickness  (t/2)  on  the  interior  walls. 

In  the  present  model  the  possibility  of  studying  the 
Trombe  type  of  wall  7  is  also  included.  Further  the  effects 
of  distributing  the  thermal  mass  between  elements  7  and  8 
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can  also  be  examined. 

The  thickness  of  the  element  8  in  the  present  study  was 
considered  to  be  5.0  inch.  This  is  based  on  the  simulation 
results  discussed  in  section  6.4. 


2 . 3  Meteoroloqica 1  Data 

Generally,  all  the  meteorological  stations  record  the 
hourly  insolation  level  on  a  horizontal  surface,  together 
with  ambient  air  temperature  and  wind  velocity  etc.  However, 
for  modelling  a  passive  house,  the  irradiation  of  vertical 
or  tilted  surfaces  is  desired.  In  section  3.1.7  an 
analytical  method  for  estimating  the  total  radiation  on  a 
tilted  surface  is  described. 

One  year's  data  "comprised  of"  the  hourly  radiation  on 
horizontal  surface  and  the  corresponding  air  temperatures 
were  used  in  the  model  as  weather  data. 


3.  THE  DIRECT  GAIN  WINDOW  MODEL 


3 . 1  Int roduct i on 

In  the  development  of  the  thermal  model  several 
possibilities  of  maximizing  the  winter  solar  gain  and 
minimizing  the  summer  overheating  were  considered.  As  a 
result,  a  more  sophisticated  method  of  calculating  the  solar 
heat  gain  through  a  single,  double  or  triple  glazing, 
employing  clear  glass  and  (or)  heat  absorbing  glass  was 
deve 1  oped . 

Calculation  of  solar  heat  gain  through  glass  has  been 
studied  by  Ye  1  lot  (16),  Parmelee,  and  et  al.  (17),  and 
Stephenson  (18).  Yel lot's  study  is  limited  to  single  glass. 
Parmelee' s  experimental  investigation  is  an  useful  source  of 
comparing  the  predicted  heat  gains  through  double  glass. 
However,  this  study  is  limited  to  double  glass. 

The  solar  heat  gain  factors  published  by  Stephenson 
(19)  and  ASHRAE  (20)  are  not  based  on  the  real  hourly 
meteorological  data.  The  equations  for  TAR  used  in  the 
calculation  of  solar  heat  gain  factors  in  Stephenson's 
method  are  based  on  a  polynomial  function.  This  method, 
although  simple,  is  not  accurate.  Another  limitation  of  the 
published  solar  heat  gain  factors  is  that  they  are 
applicable  to  a  single  pane  of  clear  glass,  with  a  fixed 
vertical  geometry. 

In  the  present  direct  gain  window  model  the  solar  heat 
gain  was  calculated  from  the  real  hourly  meteorological 
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data.  The  model  estimates  the  solar  heat  gain  through 
single,  double  or  triple  glazing  windows  employing  clear 
glass  or  heat  absorbing  glass.  It  should  be  recognized  that 
the  TAR  of  clear  and  heat  absorbing  glass  varies  with  angle 
of  incidence  of  the  beam  radiation.  Therefore,  the  TAR 
values  of  glass  were  calculated  at  each  hour  angle.  Further 
the  fraction  of  diffuse  radiation  transmitted,  and  absorbed 
is  not  the  same  as  that  of  the  direct  radiation.  Therefore  a 
method  of  calculating  TAR  of  the  glass  for  diffuse  radiation 
was  incorporated  in  the  model. 

It  is  believed  that  due  to  the  above  refinements,  the 
estimated  solar  heat  gain  will  be  relatively  accurate. 

The  solar  heat  gain  is  defined  as  the  total 
instantaneous  radiation  transmitted  through  the  glass  plus 
the  fraction  of  heat  transferred  due  to  radiation  absorbed 
by  the  glass. 

From  the  definition,  it  follows  that  the  TAR  of  a 
particular  glass  and  the  instantaneous  radiation  incident  on 
the  glass  should  be  Known  before  estimating  the  solar  heat 
gain.  This  is  illustrated  in  the  following  sections. 


3 . 2  Calculation  of  TAR  for  a  Single  Glazing 

Transmittance  of  glass  is  a  function  of  wavelength, 
angle  of  incidence  of  the  incoming  radiation,  the  refractive 
index,  and  the  extinction  coefficient  of  the  material. 

Hsieh  and  Su  (21)  have  presented  a  method  of 
calculating  the  thermal  radiative  properties  of  glass  for 
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wavelengths  between  0.3  to  206  micron.  Hsieh  and  Su' s  study 
accounts  for  the  effects  of  longwave  radiation  and  as  well 
as  the  surface  temperature  of  the  glass.  Their  results  show 
that  for  a  clear  glass  the  refractive  index  and  so  its 
extinction  coefficient  is  wavelength  independent  in 
shortwave  region  but  is  irregular  in  the  far  infrared 
region.  However,  for  most  solar  energy  applications  only  the 
shortwave  radiation  is  considered  (22).  Therefore  the  study 
of  Parmelee  (23),  which  neglects  the  effects  of  longwave 
radiation,  temperature  and  dust  etc.,  was  employed  in  the 
present  model . 

Figure  3.1  shows  a  glass  of  thickness  L  irradiated  by  a 
ray  of  intensity  I.  The  multiple  reflections  occuring  as  a 
result  of  I  are  shown  in  the  figure  3.1.,  where  '  r'  is  the 
fraction  reflected  and  'a'  is  the  fraction  of  radiation 
avialable  after  absorption  in  the  glass.  The  transmissivity 
T  is  calculated  as  follows. 

The  sum  of  the  intensities  of  the  transmitted 

components  is 


Ia(1-r)2+Ia3r2(1-r)2+Ia5r4(1-r)2+.  .  . 

This  is  a  geometric  series  and  hence  by  definition  the 
transmissivity  can  be  written  as 


T=  la ( 1 -r ) 2 [ 1+a2r 2+a4r4+ .  .  .  .  ]/I 


21 


Figure  3 . 1 


Multiple  reflections  of  direct  solar  radiation 
from  a  single  sheet  of  glass 
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on  simplification 

r  =  a ( 1  - r ) 2 / ( 1  - a 2  r 2 )  3.2. 1 

Similarly  the  reflectivity  /is  the  sum  of  the  intensities 
of  the  reflected  components  divided  by  the  intensity  of  the 
incident  radiation. 

?  =[  Ir+Ia2r ( 1 -r ) 2+Ia4r 3 ( 1 -r ) 2+ .  .  .]/! 


On  simplifying  the  sum  of  the  geometric  series  the 
reflectivity  /  of  the  single  glass  can  be  written  as 

f  =  r  +  a2r ( 1  - r ) 2 / ( 1 -a2r 2 )  3.2.2 

From  the  identity 

r+  oc  +  /  =i 

C<  =  1  -  T  -  f 

Therefore  the  absorbtivity  of  the  single  glass  can  be 
wr i t  ten  as 


c<  =  1 -r-a ( 1 -r ) 2/ ( 1 -ar )  3.2.3 

In  the  above  equations  both  '  r'  and  'a'  are  functions  of  the 
angle  of  incidence  of  the  incoming  solar  radiation. 
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The  absorption  of  radiation  in  a  partially  transparent 
medium  is  described  by  Bouger' s  law,  which  is  based  on  the 
assumption  that  the  absorbed  radiation  is  proportional  to 
the  local  intensity  in  the  medium  and  the  distance  the 
radiation  travels  in  the  medium.  By  dividing  the  sheet  of 
the  glass  into  a  number  of  layers  of  equal  thickness  dL, 
each  of  which  reduces  the  intensity  by  dl  ,  that  is 

-dl  =  I  K  dL 

Integrating  between  limits  0  to  L 

i 

(dl/I)  =  -K  L 

Io 

Log£ ( I/Io)  =  -KL 

I  -  KL 

—  =  e  =  a 

Io 


Where  'a'  is  the  fraction  of  the  energy  available  after 
passage  through  the  glass. 

The  thickness  of  the  glass  L  should  be  modified  to 
account  for  the  refraction  or  the  bending  of  the  ray.  The 
modified  L'  is  therefore  given  by 

CosOi  =  L/L' 

L '  =  L  SecOi  or  L/CosOi 
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Writing  in  terms  of  the  incident  angle 


a 


-  KL _ 

Si n20i  /  n2 


3.2.4 


The  other  term  '  r'  in  equation  3.2.1  can  be  obtained  from 
the  Fresnel's  equation  for  the  reflection  of  non  polarized 
radiation.  This  consists  of  two  components  of  radiation 
namely,  the  perpendi cu 1 ar  and  the  parallel  component  and  is 
given  by, 


Sin2  ( 0 i  -  Qr )  Tan2  (0i  -  6r ) 

2r  = - + - 

Sin2  (Oi  +  Or )  Tan2  (0i  +  Or ) 


3.2.5 


Where  the  parallel  component  is  given  by 


Sin2  ( 0i  -  Or ) 

r  =  - - 

Sin2  ( 0 i  +  Or) 


3.2.6 


The  perpendicular  component  is  given  by 


Tan2  ( 0 i  -  0r) 


r 


Tan2  ( 0 i  +  Or) 


3.2.7 
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It  should  be  recognized  that  the  above  TAR  equations 
for  single  sheet  of  glass  are  strictly  applicable  to 
monochromatic  radiation  only.  For  solar  radiation 
applications  the  average  or  total  values  of  TAR  applicable 
to  entire  solar  spectrum  have  to  be  determined. 

Deitz's  study  reported  in  reference  (22)  deals  with 
variation  of  transmissivity  of  glass  with  spectral 
wavelength.  Deitz  shows  that  the  glass  containing  little 
(0.15%  ferrous  oxide)  has  essentially  constant 
transmissivity  over  entire  solar  spectrum.  However,  the 
transmissivity  of  heat  absorbing  glass  (0.5%  ferrous  oxide) 
varies  considerably  with  wavelength. 

Parmelee  et  al .  (17)  have  made  extensive  experimental 
studies  on  the  transmission  of  solar  radiation  through 
window  glass.  For  single  sheet  of  glass  (less  than  0.15% 
ferrous  oxide)  they  found  generally  good  agreement  between 
measured  and  calculated  values. 

The  studies  of  Dietz  and  Parmelee  suggest  that  the 
transmissivity  of  window  glass  (also  called  clear  glass)  is 
constant  over  entire  solar  spectrum.  Therefore  this 
assumption  was  made  in  the  present  study. 

For  a  typical  1/4  inch  thick  clear  glass  the 
trasmi ssi vi ty  is  of  the  order  of  0.88,  the  absorptivity  of 
the  range  of  0.04,  the  rest  0.08  being  the  reflectivity. 
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3 . 3  Equa t i ons  for  TAR  for  Triple  Glazing  W i ndow  due  to 
D i rect  Radiation 

In  this  section  a  method  of  calculating  the  TAR  of 
triple  glazed  windows  is  presented.  The  general  recurrence 
formulae  developed  for  a  triple  glazed  window  can  be  applied 
to  single  as  well  as  double  glazed  windows.  It  is  believed 
that  by  applying  the  same  principles  a  general  recurrence 
relation  for  a  series  of  7 n'  glass  pane  windows  can  be 
developed.  This  however  was  not  attempted  in  the  present 
study. 

Equations  for  TAR  for  double  glazed  windows  are 
commonly  found  in  the  literature,  e.g.  see  Parmelee,  (23) 
Mitalas  &  Stephenson  (24).  Equation  of  transmissivity  for 
' n'  glass  covers  in  the  solar  collectors  has  been  reported 
(22).  However,  the  latter  case  neglects  the  absorption  of 
the  glass  and  assumes  that  the  TAR  of  each  glazing  is  the 
same . 

Recently,  equations  for  transmittance  of  '  n'  glass 
covers,  based  on  net  radiation  method  has  been  reported  by 
Speigel  (25).  Wijeysundera  (26)  has  extended  this  method  for 
transmissivity  as  well  as  absortivity  of  '  n'  glass  covers. 
However,  at  higher  incident  angles  the  results  obtained  from 
net  radiation  method  are  about  4%  higher  than  that  obtained 
by  ray  accounting  technique  (26). 

The  present  model  for  TAR  is  based  on  the  ray 
accounting  technique.  Therefore  it  is  expected  to  give  more 
accurate  results.  The  model  accounts  for  the  absorption  of 
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the  glass  as  well  as  for  the  different  TAR  values  for  each 
glazing.  There  is  therefore  enough  scope  in  that  by  changing 
the  TAR  of  each  glazing  an  effective  combination  of  clear 
glazings  can  be  designed. 

Figure  3.2  shows  a  ray  of  direct  beam  radiation  falling 
on  the  first  glazing.  Thereafter  the  multiple  reflections 
follow.  By  placing  the  second  glass  in  the  centre  the 
pattern  of  multiple  reflections  was  studied.  Let 


Ti  . 

c<  1  f 

/l 

represent 

cK  2.  i 

A 

represent 

t3, 

0^3  i 

A 

represent 

Referring  to  figure  3.2 


the  TAR  of 
the  TAR  of 
the  TAR  of 


the  first  glazing 
the  second  glazing 
the  third  glazing. 


Let 


R,  =  /iti 

Ti  =  T, 


3.3.1 

3.3.2 


It  is  clear  from  the  multiple  reflection  pattern  that 


R2  -  fi/i(  Ri  )  +  A  TT  (  Ti  ) 
T2=  /i  T2(  Ri  )  +  A/3  (Ti  ) 


Simi  1  ar  ly. 


. 
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Figure  3.2  Multiple  reflections  of  direct  solar  radiation 

from  a  triple  glazing  window. 
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R  3  ~  ^2  (  R  z  )  +  T  z  ) 

T3  =  -^Ta(Rz)  +  SzJsil-z) 

and  so  on.  Writing  R  and  T  in  general  form, 


Rj 


op 

W?£(Rj-i  )  +  4  r2(  T  j  -1  ) 

i*=  2. 

00 


Tj  = 


f,r2(  Rj-i  ) 


JU(Tj-i  ) 


T=Z 


Combining  all  the  reflected  components, 


3.3.3 

3.3.4 


/l23  =  A  + 


+  /s'G(Tj-i) 


3.3.5 


Similarly,  combining  all  the  transmitted  components, 


'"Ci  23=  'Cj  Ti  T3  +  'O. 


oO 

V  4ti(Rj-i  ) 

j=2 


A/3  (  T  j  -1  ) 


3.3.6 


From  equations  3.3.5  and  3.3.6  we  can  determine  the 
absorbtivity  o<l2 3of  a  triple  glazed  window. 


C<,Z3  =  1  -  /,23  -  Ti 2.3 


3.3.7 


cK  123  = 


oO 

f,  1-  Arf*  T,r»T,+  r,  Wrj.,1  +  /»TS(TJ-1 )) 

:r=z. 


00 

+  r, \(/<r2(Rj-i  )  +  /2£(Tj-i  D 


3.3.8 
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However  to  determine  the  net  heat  transfer  through  the 
window  the  absorbt  i  vi  t  i  es  of  the  individual  g  1  azi  ngs, 
o<Zof3  c<30f3,  should  be  Known.  The  procedure  of 
calculating  the  individual  absorbt i vi t i es  of  each  glazing  as 
an  integral  part  of  a  triple  glazing  window  is  presented 
below. 

To  determine  the  absorbtivity  of  first  glazing  of  a 
triple  glazing  window  the  components  of  the  radiation 
incident  on  the  first  glass  should  be  multiplied  by  its 
absorbt i vi ty . 

Therefore 


4r2  (Tj-i  ) 


3.3.9 


+ 


Similarly  absorbtivity  of  the  inner  glazing  is  given  by 


+ 


£>3  (  T  j  — i  ) 


3.3.10 


For  determinig  the  absorbtivity  of  the  middle  glass 


Let 
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Ui  be  the  fraction  of  the  transmitted 

component  from  the  first  glass  incident 
on  the  middle  glass. 

Bi  be  the  fraction  of  the  reflected 

component  from  the  third  glass  incident 
on  the  middle  glass. 

That  is: 


Ui  =  fi&Ti  3.3.11 

Bi  =  ?3  T1T2.  3.3.12 


From  figure  3.2  it  follows  that 


Uz  = 

Bz  = 


f,A(  Ui  )  +  fiTUBi  ) 

/3tk(Ul  )  +  /z&(Bi  ) 


Simi  1  ar  ly. 


U3 

B3 


f,A(Uz)  +  fiTz(Bz) 

F3Vo{  U  2 )  +  (  B2  ) 


and  so  on.  Writing  U  and  B  in  general  form, 


oO 


Uj  =  Wj-i)  +  fi^(Bj-i)  3.3.13 


y-2. 


CO 


Bj  = 


Y^T2(  u  j  -1  )  +  &  (  B  j  -1  )  3.3.14 


J-2 
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The  sum  of  the  intensities  of  all  components  incident  on  the 
middle  glass  multiplied  by  its  absorbtivity  gives 


°<20F3  =  (1  -  A  -  T*  ) 


'"0  +  ?\  ftz Ti  +  fjT| 'Xt. 


oO 


[UM  (  fl  +  -^3X2.)  +  Bj-1  (  ^2^3  +  f'T2 


T=2. 


3.3.15 


3 . 4  App 1 i cabi 1 i ty  of  TAR  Equat ions  of  Triple  Glazing  Window 
to  Single  and  Doub 1 e  Glazing  W i ndows 

The  transmissivity,  reflectivity,  and  absorbtivity 
equations  developed  in  the  previous  section  for  a  triple 
glazing  window  are  particularly  convenient  for  use  on  the 
digital  computers. 

The  equations  are  general  in  nature  in  the  sense  that 
they  can  be  used  for  calculating  the  TAR  of  a  single  and 
(or)  double  glazing  windows. 

The  following  analysis  demonstrate  the  applicability  of 
the  TAR  equations  for  a  double  glazing  window.  The 
transmissivity  equation  for  a  double  glazing  window  can  be 
obtained  by  substituting  T3  =  1  and  T3  =  0  in  equation  3.3.6 

cO 

T'12.  =  ~^T*+  )  +  /z-o(Tj-1  ) 

T»  2. 

oO 

-  ^r2+  1-  ( R  j-i  ) 

JrZ 


Summing  up  from  j  =  2  to  j  =  3 


33 


~C,2  =  T,T2+  ftfzTITz  (1  +  f,f 2  + 


=  t;t2 


1  +  .P,  A 


z  + 


/?  fz 


This  is  the  sum  of  an  infinite  series.  Hence  can  be  written 
as 


Ti2  =  T,Ta/  (  1  -  ) 


3.4.1 


Equation  3.4.1  is  the  equation  for  transmissivity  of  a 
double  glazing  window  (23). 

Similarly  as  an  example  the  equation  3.3.9  can  be 
reduced  to  determine  the  absorbtivity  of  the  first  glass  of  a 
double  glazing  window. 

Substituting  =  1  and  A  =  0  in  equation  3.3.9 


OC  OF  2  =  (1  -  fl  -  T5|  ) 


1  +  JzTi  + 


cO 


J-Z 


Rj 


Considering  three  terms  of  the  series  and  substituting  Ri ,  Rz, 
and  R3 


=  (1  -  /,  -  Ti  )  [  1  +  f  2  Ti  +  ft  fz.  (  -fa  Ti  )  +  -  *  • 

.  ( 1  .  ft  .  )  [  1  +  A  Ti  (  1  +  A  +/?*/£+•••)] 
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oOofz  =  (  i  -  f,  -Ti  )  [  1  -  A  (  A  -  A  )  ]  /  (1  -  fifiz  )  3.4.2 

This  is  the  equation  for  absorbtivity  of  the  first  glazing 
of  a  double  glazed  window  (23). 

Hence  in  general  it  can  be  shown  that  the  recurrence 
equations  for  TAR  of  a  triple  glazing  window  can  be  used  for 
double  as  well  as  single  glazing  windows. 

To  verify  the  accuracy  of  the  recurrence  formulae  for 
different  j  values  a  simple  arithmetic  substitution  was 
carried  out.  By  considering  the  first  two  terms  in  the 
series  the  error  in  the  calculated  value  of  transmissivity 
was  observed  to  be  1.3%.  Addition  of  another  term  improved 
the  accuracy  to  within  0.3%.  Based  on  these  observations  a 
value  of  j  =  6  was  considered  to  be  adequate.  The  TAR  values 
obtained  from  the  recurrence  equations  are  accurate  to  third 
decimal  place  with  those  reported  in  the  references  (23,24). 


3 . 5  Calculation  of  TAR  for  the  Heat  Absorbing  G 1  ass  due  to 
D i rect  Radiation 

Unlike  clear  glass  (less  than  0.15%  ferrous  oxide),  the 
heat  absorbing  glass  may  have  0.5%  or  more  of  ferrous  oxide 
which  accounts  for  most  of  absorption  in  the  glass.  The 
study  of  Dietz  reported  in  reference  (22)  shows  that  the 
transmissivity  of  heat  absorbing  glass  varies  considerably 
with  wavelength.  Hence  its  effects  cannot  be  neglected. 

Therefore  the  transmittance  of  heat  absorbing  glass  is 
a  function  of  the  spectral  wavelength,  the  angle  of 
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incidence,  and  the  KL  product  of  the  glass.  Neglecting  the 
effects  of  longwave  radiation,  temperature,  dust  and  dirt, 
the  functional  dependence  of  transmissivity  can  be  written 
as , 


T  =  f  (  A  ,  6,  K,  L  )  3.5.1 

Mi  talas  (24)  suggested  a  method  of  calculating  the 
effective  transmittance  of  window  glass  by  dividing  the 
solar  spectrum  (at  air  mass  2)  into  ten  bands  each 
containing  10%  of  solar  energy.  The  present  model  is  based 
on  this  technique  .  The  advantage  of  the  present  model  is 
that  it  is  possible  to  combine  a  single  heat  absorbing  glass 
in  a  single  or  double  glazed  windows  for  studying  the 
effective  window  designs. 

Generally  the  manufacturers  provide  the  monochromatic 
spectral  transmissivity  of  the  glass  at  normal  incidence. 
This  serves  as  an  input  to  the  present  model. 

The  technique  of  calculating  the  transmissivity  of  heat 
absorbing  glass  is  one  of  transferring  the  characteristics 
from  wavelength  dependence  to  KL  product  dependence,  and 
then  integrating  these  values  over  the  entire  spectrum.  The 
calculations  involved  are  illustrated  further. 

The  transmittance  of  glass  at  normal  incidence  was 
calculated  for  KL  values  ranging  from  0.01  to  2.0  using 
equation  3.3.6.  A  polynomial  function  of  third  degree  was 
found  to  approximate  the  above  results. 
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Tn  =  0. 9047-0. 8747 ( KL ) +0 . 3309 ( KL ) 2 -0 . 05 1 2 ( KL ) 3 

3.5.2 

A  consideration  into  the  input  (monochromatic 
transmissivity)  to  the  model  and  the  polynomial  equation 
3.5.2  suggests  that  by  dividing  the  solar  spectrum  into 
equal  energy  bands,  the  band  transmissivity  at  normal 
incidence  can  be  obtained.  That  is  at  each  band  wave 
lengths,  the  transmissivity  of  the  glass  was  obtained  by 
fitting  a  spline  function  to  the  input  data. 

From  the  band  transmissivities  it  was  possible  to 
calculate  the  corresponding  KL  values  from  the  polynomial 
equation  3.5.2.  However,  these  band  transmissivities 
correspond  to  normal  incidence.  From  the  equation  for  TAR 
given  in  section  3.3,  the  band  transmissivity  at  any  desired 
angle  of  incidence  can  be  determined.  Therefore  a  simple 
average  of  the  band  transmissivities  gives  the  spectral  or 
effective  transmissivity. 

The  transmissivity  of  the  heat  absorbing  glass  obtained 
from  the  present  model  is  accurate  to  second  decimal  point 
with  that  reported  in  reference  (24).  Typically  the 
transmissivity  of  heat  absorbing  glass  will  be  of  the  order 
of  0.20  to  0.40  and  its  absorptivity  of  the  order  of  0.60  to 


0.80. 
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3 . 6  Ca 1 cu 1  a t i on  of  TAR  due  to  D i f f use  Rad i a t i on 

In  the  preceding  sections  it  was  shown  that  the 
transmissivity  of  the  glass  is  a  function  of  angle  of 
incidence  of  direct  radiation  and  the  KL  product  of  the 
glass.  It  should  be  recognized  that  the  sky  diffuse 
radiation  comes  from  all  parts  of  the  sky. 

The  transmittance  of  diffuse  radiation  is  studied  by 
Parmelee  (23),  Mitalas  and  Stephenson  (24),  and  Puri  (27). 
Puri's  model  accounts  for  isotropy  as  well  as  anisotropy  of 
diffuse  radiation.  Anisotropy  of  the  diffuse  radiation 
introduces  further  problems.  However  by  assuming  the  sky 
diffuse  radiation  to  be  isotropic,  Parmelee  (23)  suggested  a 
graphical  integration  technique  to  determine  the 
transmissivity  of  the  glass  due  to  diffuse  radiation.  The 
present  study  extends  the  applicability  of  the  Parmelee' s 
technique  for  up  to  triple  glazing. 

From  figure  3.3  the  energy  received  by  a  vertical 
window  dq  is  given  by 

dq  =  Id  dAi  Cosi  dw 

The  solid  angle  is  given  by 

dw  =  dA2/  r2 


From  the  geometry  of  figure  3.3  dAz  can  be  written  as 


: :  5  h 
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Figure  3.3  Geometry  of  diffuse  radiation  on  a  vertical 

w i ndow . 
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dA*  =  n ( r  Sin  i )  ( r  di ) 

The  unit  energy  received  by  the  window  then  becomes 
dq/dAi  =  dEs  =  Id  Cosi  a  (r  Sini )  (rdi)/r2 
dEs  =  a  Id  Sin  ( 2 i )  di 


The  fraction  of  the  radiation  reflected  r  is  the  sum  of 
the  energy  reflected  at  each  angle  divided  by  the  total 
incident  radiation. 

Therefore  integrating  over  quarter  of  the  sky 

■njz 

0.5a  Id  r  ( i  )  S  i  n  (  2  i  )  di 


rs 


o 


0 . 5  a  Id 


rs 


n/2 

r  (  i  )  S  i  n  (  2  i  )  d  i 

SmA 

O 


3.6.1 


Similarly  the  fraction  absorbed  due  to  sky  diffuse  radiation 
'as'  can  be  written  as 


n/ 2 

''-KL//1-sin2i/n2 

as  =  e  [ Sin( 2i )  di ]  3.6.2 

V/ 

O 

Two  separate  transmissivities  have  to  be  computed  (one 
for  parallel  and  another  for  perpendicular  component  of 
radiation)  using  the  TAR  equations  derived  in  section  3.3.  A 
simple  average  will  give  the  effective  transmissivity  due  to 
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sky  diffuse  radiation.  The  same  procedure  is  applicable  for 
calculating  the  absorptivity  due  to  sky  diffuse  radiation. 

The  diffuse  transmissivity  computed  by  the  present 
model  is  accurate  to  second  decimal  place  with  that  reported 
in  reference  (23).  However,  the  predicted  values  are 
primarily  based  on  the  assumption  that  the  sky  diffuse 
radiation  is  isotropic.  The  errors  introduced  due  to  such  an 
assumption  may  not  be  large. 


3 . 7  Total  Radiation  on  a  T i 1  ted  Surface 

The  present  model  utilizes  the  hourly  insolation  on  the 
horizontal  surfaces.  But  in  the  thermal  model  the 
irradiation  of  vertical  and  (or)  tilted  surfaces  is 
required.  Several  authors  have  developed  analytical  models 
to  estimate  the  total  radiation  on  a  tilted  surface. 

Liu  and  Jordan's  model  (28)  is  perhaps  the  simplest.  It 
assumes  that  the  sky  diffuse  radiation  is  isotropic  and 
introduces  a  view  factor  to  account  for  the  tilted  surfaces. 
Parmelee  (29),  and  Threlkeld' s  (30)  experimental 
investigation  suggests  the  non- i sotropi c  nature  of  the 
diffuse  radiation.  Puri  (31)  has  developed  a  model 
considering  the  diffuse  radiation  to  be  partly  isotropic  and 
partly  non- i sot ropi c . 

Several  other  reseachers  have  studied  the  relationship 
between  the  diffuse  and  total  radiation.  Among  them  Ruth  and 
Chant  (32),  Pereira  and  Rabl  (33),  and  Hay  (34)  have 
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proposed  corrections  to  Liu  and  Jordan's  model.  Wesely  and 
Lipschutz  (35)  have  reported  a  method  of  estimating  hourly 
average  diffuse  and  direct  radiation  under  scattered  clouds. 
Nicholas  and  Child  (36)  have  published  the  solar  radiation 
charts.  All  the  above  studies  are  primarily  concerned  with 
shortwave  radiation.  However,  Cole  (37)  suggests  the  errors 
involved  in  neglecting  the  longwave  radiation  effects.  The 
study  of  Dave  (38)  shows  the  errors  involved  in  isotropic 
distribution  approximations. 

Therefore  estimation  of  total  radiation  on  tilted 
surfaces  is  a  matter  of  continuous  research.  From  the 
literature  survey,  a  combination  of  the  above  methods  was 
adopted  for  the  present  study. 

The  magnitude  of  diffuse  radiation  estimated  by  Liu  and 
Jordan's  method  for  a  tilted  surface  is  somewhat 
underestimated  (34).  It  was  observed  that  for  a  vertical 
surface,  the  HdT  calculated  by  Liu  and  Jordan's  method  was 
12-15%  less  than  that  obtained  by  ThrelKeld' s  method.  The 
difference  was  consistent.  Due  to  this  fact  the  diffuse 
radiation  on  a  tilted  surface  calculated  by  Liu  and  Jordan's 
method  was  multiplied  by  a  factor  of  1.12. 

There  is  no  justification  to  this  approach,  but  it  is 
believed  that  introduction  of  factor  like  1.12  is  never 
likely  to  overestimate  the  diffuse  radiation  component. 

Per i ved  Solar  Angles 

For  a  particular  surface  orientation,  the  sun's  incidence 
angle  0,  altitude  ,  azimuth  angle  SAZ  ,  can  be  derived 


* 
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from  three  basic  angles:  namely,  latitude(  0  ),  hour  angle 
(  w  )  and  declination  (  5  ). 

s  =  23.45  Sin  [360(284  +  n  /  365)1  3.7.8 

where  n  is  the  day  of  the  year. 

The  hour  angle  w  is  given  by 

w  =  (12.0  -  Sn)  x  15.0  degrees  3.7.2 


where 


Sn  =  ET  -1/15  (Ls  -Lg)  +  12.0 

The  Sun's  altitude  J3  ,  solar  azimuth  SAZ,  effective  azimuth 
EAZ  and  the  incident  angle  0  were  calculated  using  the 
following  equations. 


Sin^  =  Cos0  Cosw  Cos  8  +  Sin0  SinS  3.7.3 
Sin  (SAZ)  =  Cos  8  Sin  w  /  Cos#  3.7.4 
EAZ  =  SAZ  ±  HOA  3.7.5 
CosO  =  Cos yS  Cos(EAZ)  Cos  S  +  Sinp  Sin  S  3.7.6 


Neglecting  longwave  radiation  the  total  radiation  on  a 
tilted  surface  is  therefore 

Hr  = 


HbT  +  HdT  +  HrT 


3.7.7 
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Hay  (34)  showed  that  the  total  radiation  H  on  a  horizontal 
surface  can  be  split  into  beam  and  diffuse  components  with 
the  following  ratios: 

H'  H  (1  -  P&  ) 


where 


£  =  0.6C  +  (  1-0  (0.25) 


and 


Hd'  /  H'  =  1.0408  -  0.7016  (H'/Ho)  -4.4329  (H'/Ho)2 

+  6.5059  (H/Ho)3  -2.2752  (H'/Ho)5  3.7.9 


Hd'  Hd  -  H 

H'  H(  1  -  ) 


3.7.10 


From  the  above  equation  Hd  on  a  horizontal  surface  was 
ca leu  1 ated . 

Therefore 


Hb  =  H  -  Hd 


3.7.11 
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The  extraterrestrial  radiation  Ho  was  calculated  by  using 
the  equation 

Ho  =  Isc  [  1  +0.033  Cos  ( 360n/ 365 ) ]  Sinp  3.7.12 

The  magnitude  of  beam  radiation  on  the  tilted  surface  was 
estimated  as  follows. 

From  figure  3.4  Rb  can  be  written  as, 

Rb  =  CosO  /  CosOz 
OR 

Rb  =  CosO  /  Sin£ 

Therefore 

HbT=  Hb  (Rb)  3.7.13 

From  Liu  and  Jordan's  (28)  model 

Hd T  =  Hd  ( 1  +  Cos  S  /  2)  3.7.14 

and 

HrT  =  Hr  11  -  Cos  S  /  2)  3.7.15 

therefore 


Hr  =  Hb  (Rb)  +  Hd  (1  +  Cos  S/2)  +  H p  (1  -  Cos  S/2) 

3.7.16 
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Figure  3.4  Radiation  on  a  horizontal  and  tilted  surface. 
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For  a  vertical  surface  Threlkeld'  s  function  for 
estimating  Hdv  was  employed 

Threlkeld  showed  that  for  CosO  >  -  0.2 


Hdv/Hd  =  0.55  +  0.437  (CosO  +  0.313  Cos20) 


otherwi se 


3.7.17 


Hdv/Hd  =0.45 


For  surfaces  which  are  tilted,  the  total  radiation  was 
estimated  using  Liu  &  Jordan's  equation  but  with  a 
correction  factor  of  1.12 


Ht  =  Hb  (Rb)  +  1.12  Hd  (1  +  Cos  S/2)  +  Wp  (1  -  Cos  S/2) 

3.7.18 


3 . 8  Equa t i ons  for  Sun  1 i t  Fraction  of  W i ndow  due  to  D i rect 
and  Pi f f use  Rad i at i on 

An  obvious,  but  important  problem  in  solar  heat  gain 
calculations,  is  determining  the  sunlit  fraction  of  a 
surface . 

A  typical  problem  arises  when  architectural  projections 
like  solid  overhangs  are  to  be  designed.  Monthly  average 
solar  radiation  on  a  shaded  surface  is  studied  by  Utzinger 
and  Klein  (39).  For  the  present  study  a  general  case  as 
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shown  in  figure  3.5  was  considered.  The  sunlit  fraction  (of 
the  window)  due  to  direct  radiation  was  calculated  at  each 
hour  angle  from  the  geometry  of  the  overhang  and  the  sun7  s 
position  in  the  sky. 

Similarly,  the  overhang  partially  shades  the  diffuse 
radiation.  Applying  the  same  analogy  another  equation  for 
sunlit  fraction  due  to  diffuse  radiation  was  developed.  The 
applicability  of  these  equations  to  a  more  familiar  geometry 
like  a  vertical  window  is  demonstrated. 

From  figure  3.5  it  follows  that, 

X  F 

Si  n  £  Si n (  £  +  S  ) 

The  sunlit  fraction  due  to  beam  radiation  is  the  ratio  of 
the  net  sunlit  area  to  the  total  area  of  the  surface. 
Therefore 


Fb 


X 

Hw  /  Sin  ( S ) 


3.8.1 


Where 


X 


F  Sinys/  Cos£  for  a  vertical  surface. 
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Figure  3.5  Shading  of  a  sloping  window  from  direct  solar 

radiation  by  a  solid  overhang 
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The  general  equation  for  Fb  is 


F  Si n;B 

Fb  =  - 

Si n (  £+  S  ) 


E  / 

-  /  (  Hw/Sin  S 

Sin  S  / 


3.8.2 


For  a  vertical,  window  S  =  90°,  and  therefore  Fb  reduces  to 

Fb  =  [(F)  Tan/*  -  E  ]  /  Hw  3.8.3 

Similarly  the  sunlit  fraction  due  to  diffuse  radiation  can 
be  writ  ten  as 


Fd 


(1  - 


J{  Lw  Sin  S  +  E  ) 


2  +  F2  -  F 


Lw  Sin  S  +  E 


3.8.4 


Where 


Lw  Sin  (S)  =  Hw 

For  a  vertical  window  and  with  no  overhang  IS  =  90°,  F  =  0l 

yiw2 


Fd  =  1  -  ( 1 


Lw 


)  =  1 


3.8.5 
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3 . 9  Solar  Heat  Gain  Ca leu  1  at  ion 

By  definition  solar  heat  gain  is  the  sum  of  the 
radiation  transmitted  and  the  inward  flow  of  heat 
transferred  due  to  radiation  absorbed  by  the  glass. 

Solar  heat  Radiation  transmitted  Inward  flow  of  the 

=  -|- 

gain  through  the  glass  absorbed  radiation 


In  general  for  a  triple  glazing  window 


FS  =  HbT  (  Tb  +  Ni  ^Cb  )  +  HdT  (  Td  +  NicXd  )  3.9.1 


Where  Ni  is  the  fraction  of  the  absorbed  radiation 
transferred  to  the  inside  of  building  . 

For  the  direct  beam  component,  we  can  write 

FSb  =  HbT  (  Tb  )  +  Nio  (  <*bo  lo  )  +  Nim  ( <^Cbm  Im  ) 

+  Ni  i  ( c<bi  I  i  ) 
3.9.2 

Rewriting  in  a  simplified  form 

FSb  =  Feb  (HbT)  3.9.3 


Where 


< 


51 


Feb  =  Tb  +  U  oCbo  /  ho  +  U  (1/ho  +1/hi)oCbm 

+  U  (1/ho  +  1 /hi  +  1/ha)<*bi  3.9.4 


And  U  is  the  overall  heat  transfer  coefficient  of  the  window 
given  by 


1 

U  = - - - 

1/ho  +  1  /hi  +  1/hz  +  1/hi  +  Lo/kgo  +  Lm/kgm  +Li/kgi 


In  which 


3.9.5 


h  is  the  overall  heat  transfer  coefficient 

L  is  the  thickness  of  the  glass 

Kg  conductivity  of  the  glass 


In  the  above  equations  the  following  characters  denote: 


1  first  air  space 

2  second  air  space 

o  outer  glass 

m  middle  glass 

i  inside  glass 

Similarly  the  solar  heat  gain  due  to  diffuse  radiation  can 


be  wr i t  ten  as 
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FSd  =  Fed  (Hd-r  ) 


3.9.6 


Where 


Fed  =  Td  +  U  °<do  /  ho  +  U  (  1/ho  +1/hi  )  c<dm 

+  U  (1/ho  +  1  /h.  +  1/h2)<xdi  3.9.7 


Total  solar  heat  gain  is  given  by 


FS  =  FSb  +  FSd 


3.9.8 


4.  THE  THERMAL  MODEL  OF  HOUSE 


4 . 1  Introduct ion 

The  problem  of  energy  balance  of  a  house  is  one  of 
balancing  the  various  heat  gains  and  heat  losses  taking 
place  in  and  around  the  house.  In  figure  4.1  the  heat  gains 
and  the  heat  losses  through  the  various  elements  of  the 
house  are  shown. 

To  study  the  transient  behaviour  of  a  house,  the 
effects  of  different  thermal  elements  of  the  house  were 
compounded.  The  problem  then  reduced  to  one  of  solving 
transient  heat  equation  expressing  an  energy  balance  in  the 
house . 

Knowledge  of  conductances  of  the  different  elements 
together  with  heat  sources  and  ambient  air  temperatures  are 
sufficient  to  solve  these  equations.  The  procedure  is 
illustrated  in  the  sections  to  follow. 


4 . 2  Heat  Loss  through  the  Exter ior  Enclosure  E lements 

Gilpin  1977  (40)  employed  the  technique  of  classifying 
the  elements  of  the  thermal  enclosure  into  two  types, 
resistive  and  capacitive.  Employing  the  same  technique  of 
classifying  the  enclosure  elements,  a  thermal  net  work  of 
energy  fluxes,  figure  4.2,  of  a  house  was  constructed.  Each 
of  these  exterior  element  see's  a  different  temperature. 
This  temperature  is  character i zed  by  the  orientation  of  the 
element  and  is  commonly  referred  to  as  SOL-AIR  temperature. 
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Figure  4.1  Shematic  of  a  house  illustrating  the  heat  gains 

and  the  heat  losses. 


ROOM  AIR  TEMPERATURE 
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Figure  4.2  A  schematic  illustration  of  the  energy  fluxes 

to  and  from  the  room  air  temperature. 
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Therefore  it  was  assumed  that  a  North  wall  would  lose  heat 
to  SAT  (North)  and  so  on. 

The  roof  and  the  ceiling  resistances  were  combined.  It 
was  assumed  that  the  ceiling  would  lose  heat  to  SAT 
calculated  for  a  horizontal  surface. 

By  employing  the  technique  of  Boi leau  and  Latta  (41), 
the  heat  loss  through  the  below  grade  elements  was 
calculated.  It  was  assumed  that  the  below  grade  elements 
would  lose  heat  to  SAT  calculated  for  a  horizontal  ground 
surface. 

For  calculating  the  heat  loss  through  the  basement 
floor  a  reasonable  assumption  was  made  that  the  basement 
floor  loses  heat  to  an  average  ground  temperature. 

The  wood  frame  type  of  wall  construction  was  treated  as 
purely  resistive  in  nature  and  the  capacity  introduced  by 
the  plaster  board  was  lumped. 

The  heat  loss  equations,  through  each  of  the  elements 
were  wr i t  ten  as 


Qwn  =  Awn/Rwn  (Tr  -  SATn) 
Qws  =  Aws/Rws  (Tr  -  SATs) 
Qwe  =  Awe/Rwe  (Tr  -  SATe) 


4.2.1 


4.2.2 


4.2.3 


Qww  =  Aww/Rww  (Tr  -  SATw) 
Qbg  =  ^Abg/Rbg  (Tr  -  SATc) 
Qc  =  Ac/Rc  (Tr  -  SATc) 


4.2.4 


4.2.5 


4.2.6 


Qf  =  Af/Rf  (Tr-  Tg) 


4.2.7 


■ 
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The  net  heat  gain  through  the  window  is  the  sum  of  the 
solar  heat  gain  and  the  heat  loss  or  gain  through  the  window 
due  to  indoor,  outdoor  temperature  difference. 

Qwi=FS  +Awi/Rwi(Tr-Ta)  4.2.8 

The  lumped  capacity  was  treated  by  equating  the  convection 
heat  loss  from  the  body  to  the  decrease  in  the  internal 
energy  of  the  body.  Thus 

hAL  (Tl  -  Tr )  =  -  c?\  (dT  /dt) 


Rear r angi ng 


dT/dt  =  -hA L/cfV  (Tl  -  Tr)  4.2.9 

The  term  cfV/hAL  is  called  the  time  constant  of  the  lumped 
capacitance.  This  was  calculated  from  the  properties  of  the 
plaster  board.  The  heat  transfer  coefficient  in  equation 
4.2.9  accounts  for  heat  transfer  due  to  convection  as  well 
as  radiation. 

The  time  derivative  was  approximated  employing  Euler's 
me  t  hod . 


At 


Tl  +  dT/dt 
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Hence 


At 


h  Ai 
fcV 


) 


Denoting  ? cV  /  hAt.  by  RC  and  At  by  dt  the  equation  for 
temperature  Tv+1  at  some  later  point  in  time  tv+1  becomes 


"V  +  1  V  ,  ,  V  V 

Tl  =  Ti  +  dt/RC  (Tl  -  Tr) 


4.2.10 


4 , 3  Dynamics  of  the  Therma 1  Mass 

The  thick  mass  walls  were  treated  with  a  distributive 
capacitance  approach.  The  effect  of  storage  was  included  in 
the  model  to  examine  the  diurnal  swings  of  temperature.  It 
should  be  realised  that  because  of  the  capacitive  nature  of 
the  thick  mass  wall  the  temperature  and  the  heat  flow  into 
the  mass  wall  are  not  in  phase.  All  temperatures  lag  behind 
the  heat  inputs. 

To  model  the  dynamics  of  the  mass  wall  a  finite 
difference  calculation  was  used. The  mass  wall  was  divided 
into  1/2  inch  thick  layers  and  the  temperatures  at  the  nodal 
points  were  recalculated  for  each  time  step.  To  ensure 
stability  and  accuracy,  a  time  step  of  3  minutes  was 
required  for  the  finite  difference  calculations.  The 
differential  equation  together  with  the  boundary  conditions 
and  the  finite  difference  solution  scheme  are  presented 


. 
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below: 

The  transient  behaviour  of  a  finite  thickness  of  slab 
such  as  the  mass  wall  in  this  case  can  be  studied  by  the  one 
dimensional  transient,  parabolic  heat  conduction  equation: 


bl 

St 

— 

r  oC  - 

4.3.1 

bt 

a>x2 

Subject  to  boundary 

condi t ions 

dT/dt  =  0 

at  x  =  L/2 

4.3.2 

T  =  Tr 

at  x  =  0 

4.3.3 

and  the  initial  condition 

T  =  Tmin 

at  t  :  0 

4.3.4 

Where  ck  =  k/ pc  is  the  thermal  diffusivity  of  the  material. 
One  half  of  the  mass  wall  was  considered  due  to  symmetry  of 
the  problem. 

Dividing  the  mass  wall  into  M  layers  of  thickness  ax  as 
shown  in  figure  4.3  an  energy  balance  for  each  of  the  nodes 


was  wr i t  ten . 


. 
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Figure  4.3  Nodal -point  arrangement  for  thermal  mass  wall. 


Figure  4.4  Energy  balance  on  TronVbe  type  storage  wall 
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This  gives  a  set  of  finite  difference  equations. 

Let  a  t  =  time  step  in  t 

ax  =  thickness  of  each  layer. 

For  first  node  m  =  1 : 

The  energy  balance  may  be  written  as 

q  =  q,x  +  Ei 

writing  in  terms  of  rate  equations 

kA  ax  dTi 

h  As  (  Tr  -  Ti  )  =  -  (  Ti  -  Ta  )  +  /?cA  —  — 

ax  2  dt 


Applying  Euler's  method  for  approximating  the  time 
derivative  and  rearranging,  the  temperature  T-,  at  a  new  time 
step  can  be  written  as 

v  2hAt  v  v  2kAt  v  v 
T,  +  -  (  Tr  -  TT  )  -  -  (  T?  -  Tx )  4.3.5 

PC  AX  pc  AX* 

For  nodes  m  =  2  to  M  -  1 : 

For  the  central  nodes,  using  the  mathematical  formulation  of 
finite  difference  equation  in  space  coordinates,  the 
equation  for  temperature  distribution  can  be  written  as 
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dTm 

dt 


(Tm-i  -  2Tm  +  Tm+i  ) 

c<  - - - - - - 


AX  2 


Forming  an  explicit  equation  for  Tm  to  move  ahead  in  time 


Kal 

PC  AX  2 


( Tm  _i 


2Tm 


Tm-vi  ) 


4.3.6 


For  1 ast  node  m  =  M : 

The  boundary  condition  dT/dx  =  0  can  be  written  in 
terms  of  difference  equation  as 


Tm+1  “  Tm-i 
(2  Ax) 


0 


Tm+i  =  T  m-i  4.3.7 

Hence  equations  4.3.5  through  4.3.7  can  be  solved 
simultaneously  to  determine  the  surface  temperature  of  the 
storage  wall. 

In  modelling  the  Trombe  type  of  storage  wall,  the  procedure 
essentially  remains  the  same.  One  of  the  boundary  condition 
now  has  to  be  modified  to  the  fact  that  the  storage  wall  is 
oriented  to  receive  the  solar  heat  gain. 
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From  figure  4.4  the  energy  balance  can  be  written  as 

FS  -  qc  -  qK  +  E  =  0  4.3.8 

Writing  in  terms  of  rate  equations 

KA  ax  dTs 

h  As  (  Ts  -  Tr  )  =  - —  (  T,  -  T*  )  -  fcA  —  —  +  FS 

ax  2  dt 


Approximating  the  time  derivative  and  rearranging 


2h  At 
;cax 


2kAt  v  2  FS  At 

-  (  fi  -  Ta  )  +  - 

PC  AXa  fcA  AX 


4.3.9 


To  study  the  dynamics  of  the  Trombe  type  of  storage  wall  the 
equation  4.3.9  together  with  equations  4.3.6  and  4.3.7  have 
to  be  solved  simultaneously  moving  ahead  in  time. 

It  can  be  seen  that  the  heat  flows,  into  and  (or)  out 
of  the  mass  wall  at  each  time  step.  This  means  that  the 
temperatures  and  the  heat  balance  calculations  are  to  be 
recalculated  at  each  time  step. 
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4 . 4  Infi Itration  Heat  Loss 

A  great  deal  of  uncertainty  exists  in  handling  the 
infiltration  heat  loss  term.  To  date,  there  is  no  reliable 
data  which  can  be  employed  with  confidence  in  the  simulation 
studies.  Konard' s  (1)  model  for  calculating  the  infiltration 
rate  is  too  complex.  For  instance  Konard  assumes  for  a 
limited  building  configurations  a  leakage  opening 
represented  by  several  equally  spaced  holes,  and  the  mass 
flow  through  the  holes  is  computed.  The  accuracy  of  this 
model  is  yet  to  be  established.  A  simple  and  yet  reasonably 
accurate  method  of  modelling  the  infiltration  heat  loss  have 
to  be  investigated. 

However,  for  the  present  study  a  simple  approach  was 
employed  in  that  a  fixed  number  of  air  changes  per  hour  were 
assumed  to  occur.  The  average  air  change  rates  measured  by 
Dale  et  al.  (42)  on  the  experimental  solar  test  modules 
built  at  The  University  of  Alberta,  were  used.  This  air 
change  rate  is  based  on  the  total  volume  of  the  house  that 
is  including  the  below  grade  volume.  The  heat  loss  due  to  an 
air  change  rate  of  R  is  given  by 

Qinf  =  /’cVR  (Tr  -  Ta)  4.4.1 


4 . 5  The  Energy  Ba 1 ance 

Thus  far  the  heat  gains,  the  heat  losses  and  the 
storages  occur ing  in  the  model  have  been  described.  It  is  a 
matter  of  interest  to  balance  the  gains  and  the  losses  along 


65 


with  the  storages.  The  parameters  of  interest  are  the  room 
temperature  and  the  yearly  auxiliary  heat  requirement  of  the 
house.  The  yearly  auxiliary  heat  is  the  sum  of  the  furnace 
heat  at  each  fractional  time  step  to  Keep  the  house  at  a 
minimum  temperature  of  65°F.  The  upper  limit  of  room 
temperature  was  set  at  75°F.  Temperatures  beyond  the  maximum 
limit  were  recorded  as  hours  of  overheating.  Figure  4.2  is  a 
schematic  illustration  of  the  energy  fluxes  to  and  from  the 
room  air  temperature. 

The  energy  balance  equation  can  be  written  as 


Qwi  +  Qint  +Qfu  +  (Qst  +  Qst  +  Q  )  -  Qwa 

-  Qc  -  Qinf  -  Qbg  -  Qf  =  0 

4.5.1 

Expanding  the  above  terms 


As  t 

Qwi  +  Qint  +  Qfu  +  -  (Tst 

Rst 


Awa 

Tr) - (Tr  -  SAT) 

Rwa 


Awi  Ac 

-  (  Tr  -  Ta)  -  (Tr  -  SATc) 

Rwi  Rc 


Ag 

-  /cVR  (  Tr  -  Ta)  -  (Tr  -  Tg)  =  0 

Rg 


4.5.2 
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Rewr i t i ng 


Ast  Awa 

Qwi  +  Qint  +  Qfu  +  (  -  Tst  +  -  SAT 

Rst  Rwa 


Ac  Ag 

+  —  SATc  +  —  Tg  ) 

Rc  Rg 

Awi 

+  (/’cVR  +  -  )  Ta  =  C  3  ( T  r  ) 

Rwi 


4. 


Let  C  3  - 


Ast  Awa 

(  -  Tst  +  - 

Rst  Rwa 


Ac  Ag 

SAT  +  —  SATc  +  — 
Rc  Rg 


and 


Awi 

+  (/cVR  +  - ) 

Rwi 


C  2  = 


Ast  Awa  Awi  Ac  Ag 

-  +  -  +  -  +  —  +  PcVR  +  —  ) 

Rst  Rwa  Rwi  Rc  Rg 


5.3 


Tg  ) 


Ta 


4.5.4 
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Rearranging  equation  4.5.5 


Tr 


Qwi  +  Qint  +  Qfu  +  C3 
C2. 


4.5.5 


The  model  handles  three  possiblities  concerning  the 
window  shutter  system.  Therefore  the  calculation  of  the 
yearly  auxiliary  heat  requirement  differs  in  each  case. 

The  solution  scheme  is  illustrated  in  the  flow  charts  4.5 
through  4.7.  A  complete  thermal  analysis  from  the  model 
would  include  month  by  month  analysis  for  one  year,  together 
with  yearly  totals.  As  an  illustration  a  sample  computer 
print  out  for  the  month  of  January  is  shown  in  Appendix  A-4. 
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Figure  4.5  Flow  chart  illustraing  the  solution  scheme 

with  no  shutters. 
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Figure  4.6  Flow  chart  illustrating  the  solution  scheme 

with  automatic  shutters. 
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Figure  4.7 


Flow  chart  illustrating  the  solution  scheme 
with  manual  shutters. 


5.  THE  COMPUTER  PROGRAM 


In  this  section  the  flow  charts  illustrati 
scheme  of  the  thermal  model  are  presented, 
as  well  as  the  subroutines  are  depicted  in 
through  5.9 


ng  the  solution 
The  main  program 
figures  5 . 1 
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Figure  5. 


c 


START 


) 


Design  Data 


Start  SHG 
Calculations 


CALL  G  PRD  I  F 


CALL  H  BMDI  F 


CALL  HTJLT 


CALL  GPRDIR 

*  ; 


CALL  SLHGF 


5 


Start 

Thermal  Analysis 


CALL  ENRGBL 


i 


CALL  H  S  DAT  A 

* 

S'  Read 

|  Radiation  Data  H,T 

T 

Write 

Write  Monthly  Results 


1  Overall  flow  chart  for  the  main  program. 
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Figure  5.3  Flow  chart  for  subroutine  HBMDIF 


75 


Figure  5.4 


Flow  chart  for  subroutine  HTILT 
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Figure  5.7  Flow  chart  for  subroutine  SOLHGF 
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Figure  5.8  Flow  chart  for  subroutine  QINTER 
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Figure  5. 


9  Flow  chart  for  subroutine  ENRGBL 


6.  SIMULATION  RESULTS  AND  DISCUSSION 


6 . 1  Int roduct i on 

The  motivation  behind  the  development  of  the  present 
thermal  model  was  to  be  able  to  predict  the  functional 
relationship  between  a  number  of  variables  that  go  into  the 
calculation  of  the  auxiliary  heat  requirement  of  a  house. 

The  direct  gain  window  model  discussed  in  section  3.0  forms 
an  important  element  of  the  present  thermal  model.  Therefore 
a  short  discussion  on  the  results  of  the  direct  gain  window 
model  are  first  presented.  This  is  followed  by  the  more 
comprehense  results  like  the  dynamics  of  the  thermal  mass  in 
the  house,  the  general  design  curves,  and  the  economics  of 
the  passive  design  features. 


6 . 2  Tr ansmi t tance  of  W i ndow  G 1  ass  as  a  Funct ion  of  Angle  of 
Incidence 

In  the  direct  gain  passive  systems,  the  design  of 
window  element  is  most  critical.  Obviously,  a  great  deal  of 
sophistication  is  necessary  in  modelling  the  direct  gain 
window  element. 

In  the  calculations  of  solar  heat  gain,  the 
transmissivity  of  glass  due  to  beam  as  well  as  diffuse 
radiation  has  to  be  determined  accurately.  Therefore  in  the 
present  model  the  transmissivity  of  the  clear  glass  was 
calculated  at  each  hour  angle.  As  an  illustration  the 
transmissivity  of  clear  glass  as  a  function  of  angle  of 
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incidence  was  plotted  in  figure  6.1.  The  transmissivity  of 
clear  glass  was  assumed  to  be  a  function  of  angle  of 
incident  beam  radiation,  the  absorption  coefficient  of  the 
glass,  and  the  thickness  of  the  glass.  The  transmissivity  of 
clear  glass  was  assumed  to  be  constant  over  the  solar 
spectrum.  Also,  the  effects  of  longwave  radiation,  dust,  and 
the  temperature  were  neglected. 

It  is  interesting  to  note  that  the  transmissivity  of 
heat  absorbing  glass  is  22%  less  than  that  of  a  triple 
glass.  The  effectiveness  of  heat  absorbing  glass  in  reducing 
the  summer  solar  gain  is  therefore  evident.  However,  this  is 
not  a  desirable  feature  for  winter  months. 
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Figure  6.1  Transmittance  of  window  glass  as  a  function  of 

angle  of  incidence. 
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6 . 3  Effect  of  Solid  Overhang  on  the  Summer  Solar  Ga i n 

The  problem  of  minimizing  the  summer  solar  gain  yet  not 
decreasing  the  winter  solar  gain  can  be  studied  by  proper 
design  of  an  overhang  on  clear  glass  window. 

In  order  to  study  the  effectiveness  of  solid  overhang 
two  cases  were  considered.  In  the  first  case  an  overhang 
with  dimensions  of  F=2.5  feet  and  E=1.0  feet  (figure  3.5) 
was  studied  on  a  double  glazing  clear  glass  window.  In  the 
second  case  the  overhang  dimensions  were  increased  to  F  =  5 . 0 
feet  and  E=2.5  feet.  The  two  designs  of  overhang  were  also 
studied  on  a  triple  glazing  clear  glass  window. 

The  monthly  average  solar  gain  during  January  and  July 
were  plotted  for  the  above  two  cases.  This  is  shown  in 
figures  6.2  and  6.3.  It  can  be  observed  from  figure  6.2  that 
the  decrease  in  the  solar  heat  gain  due  to  a  triple  glazing 
window  is  relatively  small.  Also,  the  overhang  does  not 
shade  the  low  winter  sun  but  works  to  a  greater  advantage  in 
shading  the  summer  sun  radiation  as  shown  in  figure  6.3. 

However,  the  satisfaction  of  these  two  entirely 
opposite  requirements  from  a  single  window  element  cannot  be 
fully  realised  in  practice.  This  problem  can  be  somewhat 
elevated  by  choosing  a  window  area  appropriate  for  the 
heating  load  of  the  house.  The  set  of  design  curves 
presented  in  section  6.7  will  assist  in  this  selection. 


SOLRR  HERT  G-R  I  N  (btuh/sq-ft) 
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Figure  6.2 


Typical  solar  heat  gain  through  a  double  and 
triple  glass  window  for  the  month  of  January. 


SOLAR  HEAT  GA I N  (btuh/sq.ft) 
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Figure  6.3  Effectiveness  of  a  solid  overhang  in  reducing 

the  summer  solar  gain 
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6 . 4  Optimal  Thickness  of  Mass  Wa 1 1 

The  effect  of  increasing  the  thickness  of  mass  wall  on 
the  auxiliary  heat  requirement  of  R42-25-25  house  was 
studied.  The  house  design  parameters  are  listed  in  table  1. 
It  can  be  observed  from  figure  6.4  that  the  ideal  thickness 
of  the  thermal  mass  for  diurnal  energy  storage  is  either  in 
the  form  of  2.5  inch  thick  surfacings  or  in  the  form  of  5 
inch  thick  mass  wall  exposed  on  both  sides  to  the  room  air. 
Thickness  of  mass  wall  beyond  5  inches  does  not  improve  the 
performance  of  the  house.  Based  on  this  result  a  5  inch 
thick  mass  wall  was  considered  in  the  present  study. 


6 . 5  Optimum  Number  of  Nodes 

The  dynamics  of  thermal  mass  wall  were  studied  on 
R42-25-25  house.  Figure  6.5  shows  that  the  accuracy  of  the 
yearly  heat  requirement  calculations  depends  on  the  number 
of  nodal  points  considered  in  modelling  the  dynamics  of  the 
masonry  wall.  It  can  be  seen  that  selection  of  4  nodes  in  a 

2.5  inch  thick  masonry  surfacing  involves  an  error  of  less 
than  1%  in  the  yearly  heat  requirement  calculations. 


Table  1.  Design  Specifications  of  Houses  Modelled. 
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THICKNESS  OF  MASS  WALL  (INCHES) 


Figure  6.4  Effect  of  increasing  the  thickness  of  the  mass 

wall  on  the  yearly  heat  requirement. 


PERCENT  ERROR 
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NUMBER  OF  NODES 


Figure  6.5  Effect  of  nodal  points  on  the  yearly  heat 

requirement  calculations. 
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This  error  is  based  on  the  heat  requirement  calculations 
obtained  using  double  the  number  of  nodes  in  the  same 
thickness  of  mass  wall.  An  increase  in  the  number  of  nodes 
beyond  5  does  not  greatly  improve  the  accuracy  of  the 
calculations.  This  result  suggests  that  the  optimum  number 
of  nodes  in  modelling  the  dynamics  of  the  mass  wall  should 
be  about  5. 

figures  6.4  and  6.5  serve  as  a  useful  guide  in 
selecting  the  thickness  of  the  masonry  wall  and  the  number 
of  nodes  in  modelling  studies. 


. 
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6 . 6  Des i qn  Curves  for  the  Houses  Model  led 

Four  different  designs  with  identical  overall 
dimensions,  figure  6.6,  but  different  insulation  levels, 
table  1  ,  were  considered  for  the  present  study.  It  is 
belived  that  the  wide  variations  in  the  insulation  levels  of 
the  houses  represent  the  majority  of  the  contemporary  house 
insulation  practices. 

As  shown  in  table  1  the  houses  were  identified  by  a 
number  e.g.  R14-7- 15  for  the  first  house.  This  designation 
indicates  the  insulation  level  of  the  house,  that  is  14 
represents  the  resistance  value  of  the  ceiling,  7  that  of 
the  walls,  and  15  that  of  below  grade  basement  walls  and 
floor.  Table  1  gives  the  details  of  the  houses  such  as  area 
of  the  walls,  ceiling,  basement  (below  grade  walls  plus 
basement  floor),  their  resistance  values,  the  air  change 
rates,  and  the  heat  loss  coefficients.  The  heat  loss 
coefficients  were  calculated  considering  a  12%  window  area 
(based  on  above  grade  floor  area).  A  simple  method  of 
calculating  the  heat  loss  coefficient  of  a  house  is  given  in 
appendix  A-2.  The  present  model  assumes  a  reasonable 
internal  load  schedu le , figure  2.3,  for  a  family  of  3.35  (two 
adults  plus  two  children). 

Design  R14-7-15  was  chosen  as  a  reference  design  to 
compare  the  performance  of  other  three  designs. 


45 
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Figure  6.6  Design  dimensions  of  the  house 
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The  auxiliary  heat  requirement  was  calculated  for  increased 
south  facing  window  area  (as  a  percent  of  above  grade  floor 
area)  and  the  heat  capacity  in  the  house.  Figure  6.7  shows 
the  effect  of  increasing  the  window  area  on  the  heat 
requirement  of  the  house. 

The  relative  heating  requirement  of  the  four  houses  is 
evident  from  figure  6.7.  The  results  indicate  the  importance 
of  the  insulation  level  of  the  houses.  Considering  for  the 
moment  no  windows  and  no  internal  loads,  the  heat  required 
for  the  R70-40-40  enclosure  is  0.22  of  the  R14-7- 15 
enclosure.  In  such  a  well  insulated  house  the  internal  load 
can  supply  much  of  the  heating  requirement  so  that  the 
heating  requirement  with  internal  load  is  0.04  of  the 
R14-7-15  design.  Further  reduction  can  be  obtained  by  using 
south  facing  windows.  Eventually,  the  actual  heating 
requirement  could  be  made  essentially  zero  with  such  a 
design . 

The  curves  shown  in  figure  6.7  were  obtained  by 
assuming  the  manual  shutter  routine  discribed  in  section 
2.2.2.  It  should  be  realised  that  an  increase  in  the  window 
area  beyond  18%  of  the  floor  area  results  in  excessive 
overheating  as  shown  in  figure  6.8. 


’ 


RELATIVE  HEATING  REQUIREMENT 
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Fiqure  6.7  Effect  of  south  facing  window  area,  the  thermal 

mass  and  the  insulation  level  on  the  yearly 
heat  requirement  of  a  house. 
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Figure  6.8  Hours  of  potential  overheating  as  a  function 

south  facing  window  area. 
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6 . 7  The  General  Design  Curves 

The  curves  shown  in  figure  6.7  give  an  indication  of 
the  impact  of  various  specific  design  decisions  on  the 
heating  requirement  of  the  house.  However,  one  would  like  to 
have  a  more  general  way  of  estimating  the  effects  of  thermal 
mass,  window  area,  and  internal  heat  gain  for  a  house  of 
arbitrary  design.  For  this  purpose  a  normalization  of  the 
results  in  figure  6.7  was  sought  that  would  make  them  more 
generally  applicable  to  any  arbitrary  design  subjected  to 
different  weather  data  and  a  different  internal  load 
schedule . 

To  provide  some  guidance  in  selection  of  the 
appropriate  normalization  parameters  it  is  instructive  to 
consider  a  very  simplified  version  of  the  thermal  model  as 
shown  in  figure  6.9 

In  this  simplified  model  the  entire  thermal  mass  is 
assumed  to  be  at  the  room  temperature,  and  the  room  receives 
three  time  varying  heat  inputs  SHG(t),  I(t),  and  the  furnace 
output  Qfu.  The  heat  balance  describing  this  model  would  be 

C  dTr  /  dt  =  I  +  SHG  +  Qfu  -  U  (Tr  -  Too)  6.7.1 

To  extract  the  appropriate  normalization  parameters  from 
this  equation  each  term  or  group  of  terms  will  be  compared 
in  some  way  with  the  heat  loss  rate  U  (Tr  -  Too). 
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Figure  6.9  A  schematic  illustration  of  the  energy 

exchanges  through  the  resistive  and  capacitive 
elements  of  the  house. 
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First,  the  total  heat  energy  that  can  be  stored  in  the 
thermal  mass  is  C  (Tmax  -  Tmin  ).  This  value  should  be 
compared  with  the  typical  rate  of  heat  loss  during  the 


heat i ng 

season  which  is 

U  HDHr  /  HHr 

Where  HDHr  (Heating  degree  hours) 


And  HHr 

r' 

=  (Tmin  -  Too  )  dt  (  Too  <  Tmin  )  6.7.2 

1 

(Heating  hours) 

=  dt(Tco<Tmin)  6.7.3 

The  ratio  HDHr  /  HHr  is  just  the  mean  temperature  difference 
(Tmin  -  T^o  )  during  the  heating  season.  The  ratio  of  the 
amount  of  heat  storage  to  the  mean  rate  of  heat  loss  is  a 


storage 

time  constant , 

T 

C  (Tmax  -  Tmin) 

U  HDHr  /HHr  6.7.4 

Where 

o 

II 

>0 

O 

< 

This  is 

the  first  important  parameter  of  a  house  design  and 

it  gives  a  quantitative  assesment  of  the  effect  of  the 
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thermal  mass  of  the  house.  Appendix  A- 3  describes  a 
numerical  example  of  calculating  the  time  constant  of  a 
house . 

The  three  heat  inputs,  I,  SHG ,  and  Qfu,  to  the  house 
should  also  be  compared  to  the  heat  loss  rate  of  the  house. 
These  three  inputs  can  be  combined  and  expressed  as  only  two 
parameters  if  some  thought  is  given  to  the  matter.  Consider 
the  internal  loads.  The  internal  loads  can  be  thought  of  as 
consisting  of  a  constant  component  (Ic)  that  occurs  during 
the  day  or  at  night,  plus  a  fluctuating  component  (If)  that 
peaks  during  the  day  or  in  the  evening.  That  is 

I  =  Ic  +  If 

During  the  heating  season  the  constant  component,  Ic 
directly  replaces  part  of  the  furnace  heat  requirement  Qfu 
and  thus  it  is  appropriate  to  combine  these  two  components. 
The  fluctuating  component  can  be  thought  of  as  acting  in  a 
very  similar  way  to  the  solar  heat  gain  SHG  during  the  day 
time  and  therefore  it  is  appropriate  to  combine  (If)  with 
SHG.  In  both  cases  it  is  the  integral  over  the  heating 
season  that  is  important.  The  two  dimensionless  parameters 
that  result  are  the 


Solar  and  Internal  Gain  (SIG)  factor 
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ISHGF  x  Window  area  +  If  x  HHr 

U  x  HDHr  6.7.5 


Where 


ISHGF  = 


SHGF  dt  (  Too  -C.  Tmin  ) 


\J 

I  Year 


and  the  Annual  Heating  Load  Reduction  (AHLR)  factor 


lyeo V 


Qfu  +  Ic  x  HHr 


U  x  HDHr 


6.7.6 


The  SIG  factor  thus  gives  a  quantitative  assesment  of  the 
annual  space  heating  available  from  solar  and  the 
fluctuating  internal  loads  for  transacting  the  heat  losses. 
The  AHLR  factor  is  the  output  parameter  which  describes  how 
much  furnace  or  auxiliary  heating  is  required  in  the  house. 
The  annual  heating  requirement  Qfu,  is  related  to  the  value 
one  would  obtain  from  a  static  calculation  neglecting  the 
solar  gain,  U  x  HDHr,  through  this  latter  factor. 

Generally,  the  time  constant  of  the  house  has  been  an 
index  to  predict  the  performance  of  the  houses.  However,  the 
new  factors,  SIG  and  AHLR  along  with  the  time  constant  of 
the  house  fully  describe  the  storage  efficiency  of  a  house 
in  transacting  the  various  heat  losses  and  the  heat  gains. 
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There  are  three  meteorological  parameters  which  appear 
in  these  expressions.  These  are  the  Integrated  Solar  Heat 
Gain  (ISHG)  factor,  the  heating  degree  hour  index  (HDHr), 
and  the  number  of  heating  hours  (HHr).  The  latter  two  are 
usually  more  or  less  Known  for  a  given  location! for  Edmonton 
HDHr  =  137,000°C-hr  and  HHr  =  8042).  The  first  term,  the 
ISHG,  is  a  new  parameter  which  must  be  calculated  from 
meteorological  data.  It  describes  the  annual  useful  solar 
heat  gain  that  may  be  obtained  per  square  foot  of  window 
area. For  a  vertical,  south  facing,  double  glazed  window  in 
Edmonton  ISHG  =  182,000  Btu/ft2. 

Employing  the  above  analysis  the  effect  of  increasing 
the  SIG  factor  on  the  AHLR  factor  was  studied.  The  results 
of  figure  6.10  show  that  the  different  levels  of  the  house 
insulation  are  properly  normalized.  Several  curves  have  been 
drawn  through  the  data  to  indicate  the  general  dependence  of 
AHLR  factor  on  the  SIG  factor  for  different  values  of  time 
constant  (TC).One  would  expect  that  for  a  house  with  large 
time  constant  the  heating  requirement  would  approach  zero, 
AHLR=0,  for  SIG=1.0.  The  data  shows  that  as  the  time 
constant  is  increased  this  value  is  approached.  For  a  time 
constant  of  32  hour  which  is  considered  to  be  the  maximum 
possible  value  obtainable  in  a  passive  house  design  the  AHLR 
is  0.55.  A  more  practical  passive  house  design  might  have  a 
time  constant  of  8  hour  which  would  suggest  a  heat  saving  of 
about  25  percent  (AHLR  =0.75). 
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Figure  6.10  The  general  design  curves. 
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For  each  time  constant  there  is  a  broad  minimum  in  the  AHLR 
as  a  function  of  SIG.  This  implies  that  there  is  a  broad 
minimum  technical  optimum  window  area  for  a  given  time 
constant.  If  the  time  constant  can  be  increased  the  window 
area  can  be  increased  as  well  with  a  resulting  reduction  in 
heating  load. 


6 . 8  The  Economics  of  the  Pass i ve  Design  Features 

The  economics  of  passive  design  features  can  be  studied 
using  life  cycle  cost  analysis.  The  procedure  essentially 
remains  the  same  as  for  the  active  systems  (43).  It  should 
be  realized  that  the  following  analysis  are  intended  to  show 
the  procedure  of  studying  the  economics  of  passive  design 
feautures.  Therefore  the  numbers  arrived  at  may  not  have 
particular  significance  but  the  ideas  behind  the 
calculations  are  important. 

1.  An  estimate  of  the  additional  costs  associated  with  the 
passive  design  features  were  calculated.  The  features 
considered  were  the  shutter  system,  the  thermal  mass 
wall  and  the  costs  associated  with  improvement  in 
R-values  of  the  walls  and  the  ceiling.  The  life  of  the 
above  features  was  considered  to  be  20  years. 

2.  Over  a  period  of  20  years,  the  anticipated  costs  due  to 
repayment  of  money  borrowed,  interest,  insurance, 
operation,  maintenance  and  additional  property  taxes 
were  estimated. 

3.  Considering  a  market  discount  rate  (ie.  rate  of  return 
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on  the  best  alternative  investment),  of  12%,  an  estimate 
on  how  much  would  have  to  be  invested  now  to  meet  the 
anticipated  costs  was  made. 

4.  The  sum  of  the  investments  gives  the  present  worth  of 
the  passive  design. 

5.  The  savings  due  to  the  passive  design  features  and  a  non 
passive  house  were  compared  on  the  present  worth  basis. 

Employing  the  above  procedure,  the  following  three 
possibilities  were  examined. 

1.  The  savings  due  to  an  improvement  in  the  R-values  of  the 
enclosure  elements. 

2.  The  savings  due  to  the  installation  of  a  shutter  system, 
and 

3.  The  savings  due  to  the  installation  of  a  thermal  mass 
wall. 

Figure  6.11  shows  the  economics  of  increasing  the 
insulation  level  of  a  house.  The  R14-7-15  house  was  chosen 
as  a  reference  design.  The  extra  cost  involved  in  increasing 
the  insulation  level  from  R14-7-15  house  to  R25- 15-25  house 
was  obtained  from  reference  (44).  Similarly  an  estimate  of 
extra  cost  for  increasing  the  insulation  level  of  R14-7- 15 
house  to  R42-25-25  and  R70-40-40  was  also  obtained  from  the 
same  source.  These  cost  estimates  are  shown  in  table  2  under 
column  "Estimate  2".  Estimates  1  and  3  were  arbitrarily 
chosen  as  lower  and  higher  estimates  with  respect  to 
estimate  2. 
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F igure  6.11 


The  economics  of  the  insulation  level  of  a 
house 


Table  2.  Approximate  cost  estimates,  the  savings,  as  a  function 
of  Insulation  level  of  the  house. 
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At  each  cost  estimate,  the  additional  investment 
towards  insulation  was  compared  with  the  cummulative  fuel 
savings  over  a  period  of  20  years.  In  the  above  calculations 


the  gas  cost  was  assumed  to  be  $1.70  per  MBTU. 

For  convenience  the  four  house  designs  of  table  1  were 
identified  as  follows. 

House  Design  Number 

R 1 4-  7-15  1 

R25-15-25  2 

R42-25-25  3 

R70-40-40  4 

The  abscissa  in  figure  6.11  represents  the  above  houses 
by  their  design  numbers.  It  can  be  seen  that  as  the  cost  of 
adding  the  insulation  goes  up  the  optimum  point  shifts 
towards  insulation  level  corresponding  to  design  3. 

The  economics  of  installing  a  shutter  system  were 
studied  on  a  R42-25-25  house  with  a  time  constant  of  17 
hour.  Table  3  gives  three  cost  estimates  for  installing  a 
shutter  system.  Estimate  1  was  obtained  from  reference  (44). 
Estimates  2  and  3  were  arbtrarily  chosen.  The  window  area  as 
a  percent  of  above  grade  floor  area  was  increased  from  6%  to 
18%.  The  cummulative  fuel  savings  were  obtained  by  summing 
the  difference  between  heating  the  R42-25-25  house  with  and 
without  the  shutter  system.  The  above  analysis  were  carried 
out  at  two  gas  costs  as  shown  in  figure  6.12. 
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•% 
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Table  3 . Assumed  cost  estimates  and  the  calculated 
as  a  function  of  window  shutter  area 
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Figure  6.12  The  economics  of  shutter  system. 


The  first  and  second  gas  costs  were  assumed  to  be  $1 .70/MBTU 
and  $4. 70/MBTU  respectively.  The  $  savings  cor respondi ng  to 
two  gas  costs  are  shown  in  table  3.  The  numbers  shown  in 
paranthesis  in  table  3  correspond  to  the  second  gas  cost, 
differentiating  it  from  the  first  gas  cost  savings  which  are 
not  shown  in  paranthesis. 

It  is  evident  from  figure  6.12  that  at  a  gas  cost  of 
$4.70  per  MBTU  and  a  shutter  cost  of  $  2.5  per  square  foot, 
net  savings  would  result. 

Employing  the  above  procedure  the  economics  of 
installing  a  shutter  system  on  other  house  designs  can  be 
investigated.  It  is  expected  that  the  shutter  system  would 
be  more  economical  on  low  time  constant  houses  (  such  as 
R25- 1 5-25  )  because  of  the  fact  that  their  storage  capacity 
will  be  comparatively  small. 

The  economics  of  adding  thermal  mass  were  investigated 
on  R42-25-25  house.  The  present  cost  of  installation  of  a  5 
inch  thick  masonry  wall  is  around  $  45  /square  foot  (44). 
Table  4  shows  two  additional  cost  estimates.  Estimate  2  and 
3  were  kept  low  enough  to  investigate  the  most  economic  cost 
of  thermal  mass  wall.  The  area  of  the  distributed  thermal 
mass  (DM)  wall  was  increased  from  0.5  times  above  grade 
floor  (F)  area , ( DM= 1 / 2 F  in  table  4)  to  twice  the  above  grade 
floor  area  (DM=2F).  In  each  case  the  time  constant  of  the 
house  was  calculated.  The  cummulative  savings  were  plotted 
as  a  function  of  time  constant  of  the  house. 
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The  analysis  were  carried  out  at  two  gas  costs  as  shown  in 
figure  6.13.  The  numbers  in  paranthesis  in  table  4 
correspond  to  the  second  gas  cost  of  $  4.70  /  1 0 6 B T U  whereas 
the  other  numbers  correspond  to  the  first  gas  cost  of 
$1 . 70/ 1 06BTU . 

It  is  clear  that  savings  due  to  the  installation  of 
thermal  mass  would  only  be  possible  at  a  cost  as  low  as  $  1 
per  square  foot.  In  such  an  event  the  optimum  time  constant 
of  the  house  and  therefore  the  area  of  mass  wall  can  be 
selected  from  figure  6.13. 

The  above  analysis  are  essentially  intended  to  show  the 
procedure  of  studying  the  economics  of  passive  house 
designs.  The  three  cost  estimates  were  selected  to  account 
for  a  wide  range  of  uncertainties.  It  is  believed  that 
somewhere  between  the  maximum  and  the  minimum  estimates  the 
actual  figures  should  lie.  Therefore  the  trends  shown  in 
figures  6.11  through  6.13  would  be  applicable  to  a  majority 
of  the  situations. 
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TIME  CONSTANT  (HOUR) 


Figure  6.13  The  economics  of  thermal  mass  wall 


7.  CONCLUSIONS  AND  RECOMMENDATIONS 


7 . 1  Cone  1  us i ons 

Some  of  the  important  conclusions  drawn  from  this  study 

are  as  fol lows : 

1.  For  diurnal  energy  storage  the  ideal  thickness  of  a 
masonry  wail  should  be  either  in  the  form  of  2.5  inch 
thick  surfacing  or  in  the  form  of  5  inch  thick  mass  wall 
exposed  on  both  sides  to  room  air. 

2.  For  stability  and  accuracy  of  the  finite  difference 
calculations  a  minimum  of  A  nodes  or  preferably  5  nodes 
have  to  be  considered  in  studying  the  dynamics  of  the 
masonry  wall.  This  suggests  that  in  a  5  inch  thick 
masonry  wall  or  alternatively  2.5  inch  thick  surfacings 
the  thickness  of  each  layer  should  be  1/2  inch. 
Accordingly  the  time  step  works  out  to  be  3  minutes. 

3.  An  important  consideration  in  the  design  of  a  passive 
house  is  sizing  the  south  facing  windows.  The  SI6  and 
the  AHLR  factors  together  with  the  general  design  curves 
given  in  figure  6.10  are  the  useful  means  of  achieving 
this  objective. 

4.  A  comparative  study  on  four  houses  (table  1  employing 
life  cycle  cost  analysis  suggests  that  the  most  economic 
insulation  level  of  a  house  would  be  comparable  to  that 
of  R42-25-25  design. 

5.  In  order  for  the  shutter  system  to  be  economical  on  a 
R42-25-25  house  (time  constant  17  hour1,  the 
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installation  cost  has  to  be  around  $2. 5/square  foot. 

6.  In  order  for  the  thermal  mass  wall  to  be  economical  a 
cost  reduction  factor  of  about  30  has  to  be  achieved 
over  its  present  installation  cost. 

7 . 2  Recommenda  t i ons 

1.  The  two  areas  requiring  further  consideration  are  the 
modelling  techniques  for  infiltration  and  the  basement 
heat  loss  calculations. 

2.  The  accuracy  of  the  predictions  made  by  the  present 
thermal  model  has  to  be  verified  by  carrying  out 
experimental  studies. 

3.  The  economic  analysis  and  the  wood  frame  construction 
practices  suggest  that  alternative  storage  methods  have 
to  be  examined  to  replace  the  heavier  and  more  expensive 
masonry  walls. One  possibility  is  to  study  the  energy 
storage  using  eutectic  mixtures. 

4.  An  extension  of  this  study  would  be  to  look  at  the 
thermal  model  from  both  heating  as  well  as  cooling  point 
of  view  using  passive  concepts. 
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APPENDIX  A- 1 
Conversion  Table 


Quant i tv 

Enq 1 i sh  units 

Length 

1 

Ft 

Area 

1 

FT2 

Vo  1 ume 

1 

FT3 

Mass 

1 

Lb 

Densi ty 

1 

Lb/Ft3 

Specific  Heat 

1 

Btu/Lb-F0 

Insolation 

1 

Btu/Sq.ft-Hr 

Heat  Transfer 

1 

Btu/Sq. ft-Hr-F0 

Coefficient 

Therma 1 

1 

F°-Hr-Sq. ft/Btu 

Resistance 

Heat 

1 

Btu 

1 

MB  tu 

Energy 

1 

Btu/Hr 

Heat  Loss 

1 

Btu/Hr-F0 

Coefficient 

Heat  Capacity 

1 

Btu/F° 

SI  Units 
0 . 3048  m 
0.0929  m2 
0.0283  m3 
0.4535  Kg 
16.018  Kg/m3 
4. 1868  KJ/Kg-C0 
0.088  K J/m2 -Hr 
5.678  W/m2 -C° 

0. 1 76 1m2 -C° /W 

1 .055  KJ 
1 .055  GJ 
0.293  Watt 
0.5274  W/C° 

1 .9  K J/C° 


<n’ 
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APPENDIX  A-2 

Calculation  of  Heat  Loss  Coefficient  (U)  of  House 


Considering  R42-25-25  house 
From  table  1 


Awa 

= 

1209 

Sq.ft 

Rwa 

- 

25 

F°-Hr-Sq. ft/Btu 

Abg 

- 

2136 

Sq.ft 

Rbg 

- 

25 

F°-Hr-Sq. ft/Btu 

Ac 

= 

1337 

Sq.ft 

Rc 

- 

42 

F°-Hr-Sq. ft/Btu 

Awi 

- 

151 

Sq.ft 

Rwi 

- 

14 

F°-Hr-Sq. ft/Btu 

R 

- 

0.25 

( di mens i on  1  ess ) 

capaci ty  of 

air  can  be  taken  as  0.018  Btu/F° 

Total  volume  of  house  (above  as  well  as  below  gradelwas 
estimated  to  be  20160  Ft3. 

Neglecting  the  ground  temperature  effects  and  assuming 

resistance  value  of  shutter  in  closed  position,  the  heat 

loss  coefficient  was  obtained  as  follows 

Awa  Abg  Ac  Awi 
=  -  +  -  +  —  +  -  +  fc  V  R 


Rwa  Rbg  Rc 


Rwi 


1209  2136  1337  151 

-  +  -  +  -  +  -  +  0.018  (20160)  (0.25) 

25  25  42  14 


267  Btu/F° -Hr 
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APPENDIX  A-3 

Calculation  of  Time  Constant  of  House 
Considering  R42-25-25  house 

From  table  1  the  heat  loss  coefficient  U  =  267  Btu/F°-Hr 
Let  the  house  has  5  inch  thick  mass  wall  of  total  surface 
area  equal  to  630  sq.ft 

Heat  capacity  of  the  mass  wall  is  equal  to 

=  fcV 

120  (0.24)  (5/12)  (630) 

=  7560  Btu/F° 

From  equation  6.7.4  the  time  constant  of  a  house  is  equal  to 

C  (Tmax  -  Tmin) 

r-  - 

U  HDHr/HHr 

For  Edmonton 
HDHr  =  10286  (24) 

HHr  =  8042 

Wi th  Tmax  -  Tmi n  =  10° 

7560  ( 10)  (8042) 

267  (10286)  (24) 

Therefore  T  =9.2  Hour 
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APPENDIX  A-4 
Sample  Computer  Run 


ENGINEERING  DEPARTMENT 
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